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EXECUTIVE  SUMMARY 

The  volatiliiation  of  organic  chemical  compounds  of  public  health  concern  from 
streams  and  rivers  has  been  investigated.  The  volatilization  rates  are 
expressed  in  terms  of  liquid  film  coefficients  (Ki_).  The  main  objectives 
of  the  study  are:  (1)  To  determine  liquid  film  coefficients  for  ethylene  gas 
using  field  data;  (2)  To  develop  relationships  among  liquid  film  coefficients, 
channel  hydraulic  parameters  and  properties  of  chemicals;  and  (3)  To  summarize 
relation  between  chemical  structure  and  volatilization  for  organic  compounds  of 
concern  in  provincial  streams  and  rivers. 

Factors  affecting  volatilization  of  organic  compounds,  and  conceptual  models  of 
transport  across  the  air-water  interface  have  been  presented.  The  liquid  film 
coefficients  of  ethylene  are  determined  from  data  on  ethylene  and  rhodamine  WT 
dye  distributions  collected  by  the  Water  Resources  Branch  in  the  following 
streams  and  rivers: 

Avon  River  below  Stratford  Water  Pollution  Control  Plant  (WPCP) 
Baden  Creek  below  Baden  WPCP 
Grand  River  near  Woolwich  Township 
Nith  River  below  Wellesley  WPCP 
Speed  River  near  Puslinch  Township 

The  range  of  channel  hydraulic  parameters  and  liquid  film  coefficients  of 
ethylene  gas  for  the  survey  conditions  in  these  streams  are  as  follows: 

Discharge:  0.0453  -  7.29  m3/s 

Channel  width:  3.11  -  38.6  m 

Average  depth:  0.092  -  1.09  m 

Mean  velocity:  0.028  -  0.27  m/s 

Channel  bed  slope:  0.069  -  0.635  m/km 

Liquid  film  coefficients  of  ethylene  gas  (Kl)  :  0.97  -  8.07  m/day 

Temperature:  16°  -  25°C 

Relationships  among  <l,  channel   hydraulic   parameters   and   chemical 
properties   were  developed  through  statistical  and  dimensional   analysis 
procedures.  The  relationships  for  K|_  applicable  to  ethylene  gas,  given  by 


Equations  4.9  and  $.10  which  are  derived  from  dimensional  analysis,  have  beer 
found  to  provide  better  correlations  with  the  observed  Kj_  values.  The 
computed  values  of  K^  for  ethylene  gas  for  the  above  provincial  streams  are 
1n  the  range  1.87  -  5.57  m/day  from  Equation  4.9,  and  2.51  -  3.61  m/day  from 
Equation  4.10.  The  parameters  included  in  Equation  4.9  are  width,  depth, 
velocity,  drop  in  height  in  a  reach,  density  and  viscosity  of  river  water,  and 
molecular  diffusivity  of  ethylene  gas,  whereas  Equation  4.10  includes  depth, 
velocity,  density,  vSscosity  and  diffusivity.  A  sensitivity  analysis  of 
various  parameters  appearing  in  Equations  4.7-4.10  has  been  carried  out. 

For  other  streams  for  which  data  are  reported  by  Rathbun  and  Ta1  (1981), 
Equation  4.10  is  found  to  overestimate  the  K[_  values  (by  a  factor  of  3.53 
on  an  average  basis).  The  statistical  relationship  developed  by  Rathbun  and 
Tai  (1982)  in  which  JCj_  is  a  function  of  depth,  velocity  and  bed  slope,  has 
been  found  to  underestimate  K[_  values  for  the  above  provincial  streams  and 
rivers  (by  an  average  factor  of  0.75).  These  descrepancies  are  believed  to  be 
due  to  the  effects  of  unaccountable  factors  such  as  wind  velocity  and  presence 
of  other  chemicals  in  the  rivers.   It  is  suggested  that  the  relationships 
derived  from  dimensional  analysis  are  suitable  to  estimate  K|_  for  ethylene 
gas  in  provincial   streams  and  rivers.   However,  modifications  of  these 
relationships  to  include  the  effects  of  wind  velocities  is  desirable,  as  data 
on  these  aspects  becosje  available. 

An  evaluation  of  the  relations  between  Kj_  and  molecular  properties  of 
volatile  compounds  has  been  carried  out  with  the  aid  of  theoretical  review 
supplemented  by  experimental  data  reported  in  the  literature.  This  evaluation 
has  revealed  that  the  relation  between  K|_  and  critical  volume  (or  molecular 
diameter)  is  suitable  to  calculate  <i   values  for  highly  volatile  organic 
chemicals  (Henry's  law  constant  >10~3  atm-m3/mol)  using  the  values 
for  ethylene  gas. 

The  computation  of  K-l  for  a  given  organic  compound  using  the  relationships 
developed   herein  Is  outlined  in  a  step-by-step  procedure.   An  example 
illustrates  the  various  computations  involved.   In  general,  the  method  is 
applicable  to  volatile  compounds  with  Henry's  law  constants  greater  than 
10" 3  atm-m3/mol. 
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The  relationship  between  volatilization  and  chemical  structure  for  a  number  of 

organic  compounds  has  been  evaluated  based  on  the  Henry's  law  constants  (H). 

Compounds  with   H   values   <10"7  atm-mJ/mol   are   nonvolatile, 

10~7   to  10" J  atm-m3/mol  are  moderately  volatile   and 

>10"J  atm-m'/mol  are  highly  volatile.   The  organic  compounds  of 

concern  in  the  province  have  been  grouped  into  a  number  of  classes  that  are 

expected  to  behave  similarly  with  respect  to  volatilization-   Methods  of 

computing  Henry's  law  constants  and  molar  and  critical  volumes  for  organic 

compounds  have  been  presented. 
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On  a  etudie  la  volatilisation  des  composes  organiques  presents 
dans  les  cours  d'eau  et  qui  posent  des  risques  pour  la  sante 
publique.   Les  taux  de  volatilisation  sont  exprimes  sous  forme  de 
coefficient  de  transfert  de  masse  (KL) .   Cette  etude  visait 
principalement  :  1)  a  determiner  les  coefficients  de  transfert  de 
masse  de  1' ethylene  a  partir  de  donnees  techniques;  2)  a  elaborer 
les  relations  entre  les  coefficients  de  transfert  de  masse,  les 
parametres  hydrodynamiques  des  canaux  (vitesse,  profondeur  et 
largeur  du  cours  d'eau)  et  les  proprietes  des  composes 
organiques;  et  3)  a  etablir  les  relations  entre  la  structure 
chimique  et  la  volatilisation  des  composes  organiques  qui 
constituent  un  danger  pour  la  sante. 

Le  rapport  presente  les  facteurs  qui  influent  sur  la 
volatilisation  des  composes  organiques  ainsi  que  la 
conceptualisation  de  leur  transport  dans  1' interface  air-eau. 
Les  coefficients  de  transfert  de  masse  de  1' ethylene  ont  ete 
determines  a  partir  de  donnees  sur  la  distribution  de  1' ethylene 
et  du  colorant  rhodamine  WT  qui  ont  ete  etablies  par  la  Direction 
des  ressources  atmospheriques  pour  les  cours  d'eau  ci-apres  : 

riviere  Avon,  en  aval  de  la  station  de  depollution  de 

l'eau  de  Stratford 

ruisseau  Baden,  en  aval  de  la  station  de  Baden 

riviere  Grand,  pres  du  canton  de  Woolwich 

riviere  Nith,  en  aval  de  la  station  de  Wellesley 

riviere  Speed,  pres  du  canton  de  Puslinch 

Les  echelles  des  trois  parametres  hydrodynamiques  (vitesse, 
profondeur  et  largeur)  et  des  coefficients  de  transfert  de  masse 
de  1' ethylene  au  moment  de  1' etude  etaient  les  suivantes  : 

Debit  :  0,0453  -  7,29  m3/s 

Largeur  du  cours  d'eau  :  3,11  -  3  8,6  m 

Profondeur  moyenne  :  0,092  -  1,09  m 

Vitesse  moyenne  :  0,028  -  0,27  m/s 

Pente  du  canal  :  0,069  -  0,635  m/km 

Coefficients  de  transfert  de  masse  de  1' ethylene  (KL)  : 

0,97  -  8,07  m/jour 

Temperature  :  16°    -  25  *C 

Grace  a  des  analyses  statistiques  et  dimensionnelles,  on  a  pu 
elaborer  les  relations  entre  le  Kl,  les  parametres 
hydrodynamiques  des  canaux  et  les  proprietes  chimiques  des 
composes  organiques.   Les  relations  entre  le  Kl  et  1' ethylene, 
derivees  des  equations  4.9  et  4.10  —  resultant  elles-memes  d'une 
analyse  dimensionnelle  — ,  donnent  de  meilleures  correlations 
avec  les  valeurs  observees  du  KL.   La  valeur  du  KL  de  1* ethylene 
dans  les  cours  d'eau  enumeres  ci-dessus  va  de  1,87  a  5,57  m/jour 
avec  1' equation  4.9  et  de  2,51  a  3,61  m/jour  pour  1' equation 
4.10.   Dans  1' equation  4.9,  on  a  utilise  les  parametres 


suivants  :  largeur,  profondeur,  vitesse,  denivellation  du 
troncon,  densite  et  viscosite  de  l'eau  ainsi  que  diffusion 
moleculaire  de  1' ethylene.   Par  contre,  dans  1* equation  4.10,  on 
n'en  a  utilise  que  cinq  :  profondeur,  vitesse,  densite,  viscosite 
et  diffusion  moleculaire.   De  plus,  une  analyse  de  sensibilite  de 
divers  parametres  figurant  dans  les  equations  4.7  a  4.10  a  ete 
effectuee. 

Dans  le  cas  des  autres  cours  d'eau,  mesures  par  Rathbun  et  Tai  en 
1981,  on  a  constate  que  1' equation  4.10  exagerait  les  valeurs  du 
Kl  (d'un  facteur  de  3,53  en  moyenne) .   Par  ailleurs,  le  rapport 
statistique  calcule  par  les  memes  auteurs  en  1982,  selon  lequel 
le  Kl  est  fonction  de  la  profondeur,  de  la  vitesse  et  de  la  pente 
du  canal,  sous-estime,  d'un  facteur  moyen  de  0,75,  les  valeurs  du 
Kl  pour  les  cours  d'eau  qui  font  1'objet  de  la  presente  etude. 
On  croit  que  ces  ecarts  sont  causes  par  l'effet  de  facteurs  non 
traites  tels  la  vitesse  du  vent  et  la  presence  d1 autres 
substances  chimiques  dans  les  cours  d'eau.   II  semble  done  que 
les  relations  derivees  de  1' analyse  dimensionnelle  pennettraient 
d'estiiner  le  Kl  de  1' ethylene  dans  les  cours  d'eau  de  la 
province.   Cependant,  il  serait  bon  de  les  modifier  en 
incorporant  l'effet  de  la  vitesse  du  vent,  au  fur  et  a  mesure  que 
l'on  recueille  les  donnees  a  cet  egard. 

Une  evaluation  des  relations  entre  le  Kl  et  les  proprietes 
moleculaires  des  composes  organiques  a  ete  effectuee  a  l'aide 
d'un  examen  theorique  et  des  resultats  d' experience  qu'on  a  pu 
relever  dans  la  documentation.   Cette  evaluation  a  revele  que  la 
relation  entre  le  Kl  et  le  volume  critique  (ou  diametre 
moleculaire)  permet  de  calculer  les  valeurs  du  Kl  des  composes 
organiques  tres  volatils  (constante  de  la  loi  de  Henry 
>10~J  atm-m3/m°l)  er>  utilisant  les  valeurs  de  1' ethylene. 

Le  calcul  du  Kl  d'un  compose  organique  donne  a  partir  des 
relations  elaborees  dans  la  presente  etude  est  decrit  etape  par 
etape.   Un  exemple  precis  illustre  les  differents  calculs 
effectues.   En  general,  il  est  possible  d'appliquer  la  methode  de 
calcul  aux  composes  organiques  dont  les  constantes  de  la  loi  de 
Henry  sont  superieures  a  10~3  ata-xa* /ttol . 

La  relation  entre  la  volatilisation  et  la  structure  chimique  d'un 
certain  nombre  de  composes  organiques  a  ete  evaluee  a  l'aide  des 
constantes  de  la  loi  de  Henry  (H) .   Des  constantes  H  de 
<10~7  atm-m3 /mol  representent  des  composes  non  volatils;   des 
constantes  H  de  10" '  a  10~3  atm-mJ /mol,  des  composes  moderement 
volatils;  et  des  constantes  de  >10~5  atm-M*/mol,  des  composes 
tres  volatils.   Les  composes  organiques  qui  presentent  un  danger 
pour  la  sante  publique  en  Ontario  ont  ete  regroupes  en  un  certain 
nombre  de  classes,  selon  la  facon  dont  ils  reagiront  probablement 
a  la  volatilisation.   Par  ailleurs,  on  presente  ici  les  methodes 
de  calcul  des  constantes  de  la  loi  de  Henry  ainsi  que  les  volumes 
molaire  et  critique  des  composes  organiques. 
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RECOMMENDATIONS 

During  this  investigation,  several  aspects  related  to  development  of  predictive 
relationships  for  volatilization  from  water  bodies  were  found  to  require 
further  studies.  Based  on  this  knowledge,  the  following  recommendations  are 
made : 


1.  Development  of  modified  relationships  for  K^_  to  include  wind  velocity 
(through  dimensional  analyses)  is  recommended. 

2.  An  evaluation  of  the  applicability  of  the  relationships  developed  herein 
for  chemical  compounds  other  than  ethylene  is  recommended  as  data  become 
available. 

3.  It  is  recommended  that  field  data  gathering  for  volatilization  studies 
include  channel  hydraulic  parameters  (discharge,  width,  depth,,  velocity 
and  bed  slope),  and  wind  velocity  measurements. 

4.  Investigation  of  relationships  between  <l  and  molecular  properties  for 
chemical  compounds   with   Henry's    law   constants    lower    than 
10"J  atm-m3/mol  (i.e.,  moderate  and  low  volatile  compounds)  is 

also  recommended. 
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I    INTRODUCTION 

1.1   Background 

Assessment  of  potential  public  health  and  other  risks  of  organic  chemicals  in 
streams  and  rivers  requires  detailed  evaluations  of  their  transport*  cycling  and 
fate  within  the  aquatic  environment.  Therefore,  considerable  importance  is  toeing 
given  to  gathering  data  on  various  physical,  chemical  and  bioTogicaT  processes  and 
factors  related  to  these  organic  chemicals  of  public  health  concern.  One  of  the 
important  processes  common  to  many  of  these  organic  compounds  fa  volatilization. 
It  is  particularly  necessary  to  know  the  volatilization  rates  of  chemical  compounds 
when  dealing  with  emergency  situations  such  as  accidental  spills  of  suchr  compounds 
1n  order  to  predict  the  transport,  fate  and  public  health  risks  of  the  spills*. 

The  volatilization  process  represents  the  physical  transport  of  organic  compounds 
through  the  water-air  interface  into  the  atmosphere.  And,  hence,  there  is  a 
similarity  between  the  widely-known  atmospheric  reaeration  and:  the  volatilization 
processes.  Recently,  Mackay  and  Yuen  (1980)  and  Rathbun  and:  Tat  (1982)  have  shown 
that  the  volatilization  rates  applicable  to  organic  compounds  can  he  estimated  from 
the  field  data  gathered  for  the  measurement  of  reaeration  coefficients  using  dye 
and  ethylene  as  tracers. 

The  Water  Resources  Branch  of  the  Ontario  Ministry  of  the  Environment  has  conducted 
reaeration  measurement  studies  in  many  streams  and  rivers  using  dye  and  ethylene. 
These  data  can  thus  be  used  to  determine  the  volatilization  rates  applicable  to 
streams  and  rivers.  The  field  measurement  technique  is  documented  rn  thev  MOE 
publication  Water  Resources  Paper  13  (A.  Bacchus,  1981). 

1.2   Volatilization  Rate  and  Liquid  Film  Coefficient 

Volatilization  rate,  also  known  as  "desorption  rate",  is  a  coefficient  whicSi  is 
analogous  to  the  first-order  rate  coefficient,  expressed-  in  the  tinits 
day"1.  The  liquid  film  (or  mass  transfer)  coefficient,  which,  is  a  more 
fundamental  quantity,  is  utilized  in  the  expression  for  flux  across  the  axr-water 
Interface  (see  Equation  2.1  in  the  next  chapter);  and  it  is  expressed  in  the  units 
m/day.  The  volatilization  rate  (Kv)  and  liquid  film  coefficient  (KjJ  are 
related  by 
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where  I*  Is  a  characteristic  length,  equal  to  the  thickness  of  liquid  film 
layer  (see  Section  2.2). 

The  thickness  of  the  liquid  film  coefficient  can  be  looked  upon  as  being  somewhat 
analogous  to  that  for  a  boundary  layer.  In  turbulent  flow  regimes  normal  y 
encountered  in  natural  streams  and  rivers,  the  boundary  layer  thickness  is  commonly 
defined  as  the  magnitude  of  the  normal  distance  from  the  boundary  surface  i.e., 
channel  bed)  at  which  the  velocity  is  equal  to  99%  of  the  surface  velocity  (Chow, 
1959)  Thus,  for  all  practical  purposes,  the  thickness  can  be  approximated  by  the 
depth'of  water  in  the  channel.  Based  on  this  analogy,  the  characteristic  length, 
L*  is  set  equal  to  the  channel  depth,  Z;  such  an  approximation  has  been 
utilized  by  Rathbun  and  Tai  (1981).  Setting  L*=Z  in  Equation  1.1,  we  get 

KL=ZKV  ^ 

In  this  investigation,  KL  values  are  utilized  in  various  analyses  because  of 
its  fundamental  nature. 

1.3   Objectives 

The  main  objectives  of  this  study  are  as  follows: 

(a)  To  determine  mass  transfer  coefficients  applicable  to  organic  compounds  of 
public  health  concern  using  the  readily  available  data  of  ethylene  and  dye 
concentration  distributions  from  provincial  streams  and  rivers. 

(b)  To  develop  relationships  between  mass. transfer  coefficients  and  bulk-flow 
characteristics  of  river  channels  which  could  be  utilized  to  assess 
volatilization  of  organic  compounds  in  other  provincial  streams  and  rivers. 

(c)  To  summarize  relationships  between  chemical  structure  and  volatilization 
for  compounds  of  potential  concern  to  the  province. 


H    THEORETICAL  BACKGROUND 
2.1   Factors  Affecting  Volatilization 

The  volatilization  of  an  organic  compound  from  the  water  phase  to  the  air  phase  is 
dependent  on  the  physical  and  chemical  properties  of  the  compound,  the  presence  of 
other  chemical  compounds,  the  hydrodynamic  and  other  physical  properties  of  the 
water  body,  and  the  physical  properties  of  the  atmosphere  above  the  water  surface 
(Liss  and  Slater,  1974;  Mackay  and  Wolkoff,  1973;  Lyman  et  al ,  1982).  These 
factors  are  presented  below. 

The  physical  and  chemical  properties  of  organic  compounds  affecting  volatilization 
Include  molecular  diameter,  molecular  weight,  Henry's  law  constant  and  diffusion 
coefficient.  The  volatilization  rates  are  also  influenced  by  the  presence  of  some 
modifying  materials  including  adsorbents,  electrolytes,  emulsions  and  organic  films 
*  (Mackay  and  Leinonen,  1975;  Cohen  et  al,  1978). 

The  hydrodynamic  and  other  physical  properties  of  the  river  channel  include  flow 
rate,  width,  depth,  velocity,  bed  slope,  bed  roughness,  turbulence  level  and 
wind-induced  currents.  Factors  such  as  suspended  sediment  concentration  and  the 
presence  of  other  chemicals  in  the  water  phase  could  also  affect  volatilization. 

The  atmospheric  properties  of  concern  include  wind  speed,  stability  and  other 
factors.  Temperature  affects  the  vapour  pressure,  and  solubility,  and  thus, 
influences  volatility  of  chemicals. 

2.2   Theoretical  Aspects 
2.2.1  The  Two-Film  Theory: 

Theoretical  concepts  of  volatilization  of  chemical  compounds  from  water  to  the 
atmosphere  have  been  presented  by  various  researchers  (Lewis  and  Whitman,  1924; 
Mackay  and  Wolkoff,  1973;  Liss  and  Slater,  1974;  Mackay  and  Leinonen,  1975;  Chiu 
and  Freed,  1977  and  1979;  Smith  and  Bomberger,  1979;  Smith  et  al,  1980).  A 
review  of  these  theoretical  approaches  and  their  limitations  have  been  presented  in 
a  recent  publication  (Lyman  et  al,  1982).  The  two-film  theory  of  Lewis  and  Whitman 
(1924)  is  generally  utilized  to  describe  the  volatilization  of  organic  compounds 


from  water.  The  atmospheric  reaeration  in  river  channels  is  also  described  by  the 
two-film  theory.   The  salient  aspects  of  the  two-film  theory  are  presented  below. 

The  two-film  model  assumes  that  the  bulk  air  and  water  phases  are  uniformly  mixed, 
and  that  these  two  phases  are  separated  by  thin  films  of  air  and  water,  as  shown 
schematically  in  Figure  2.1.  The  main  resistance  to  gas  transport  is  considered  to 
exist  in  the  liquid  and  gas  phase  interfacial  layers  (or  films).  Within  these  two 
films,  the  transport  of  the  exchanging  gas  is  assumed  to  take  place  by  molecular 
diffusion.  Then,  application  of  the  Fick's  first  law  of  diffusion  for  gas 
transport  through  each  layer  results  in 

F  =  kL  (CSL  -  CJ  -  kg  (CG  -  Csg)  (2.1) 

where  F  is  the  flux  of  gas;  k|_  and  kQ  are  the  mass  transfer  coefficients 
for  the  liquid  and  gas  phases,  respectively  (meters/day):  C|_  and  Cq  are 
the  concentrations  in  the  bulk  liquid  and  gas  phases,  respectively;  and  C$[, 
and  C$g  are  ^ne  liquid-phase  and  gas-phase  concentrations  at  the  interface, 
respectively  (see  Figure  2.1).   The  mass  transfer  coefficients  (also  termed 
"exchange  constants")  are  defined  by  the  general  relation 

k   =  Dm/6  (2.2) 

where  0m  is  the  coefficient  of  molecular  diffusion  of  gas  in  the  appropriate 
layer  of  thickness  6  (see  Figure  2.1).  The  coefficients,  k|_  and  kg, 
are  measures  of  the  flux  of  gas  per  unit  concentration  gradient. 

Following  Cohen,  et  al  (1978),  the  transport  across  the  two-layer  system  (shown  in 
Figure  2.1)  can  also  be  expressed  as  follows  on  the  assumption  that  the 
concentrations  immediately  on  either  side  of  the  interface  are  in  equilibrium  as 
expressed  by  a  Henry's  law  constant  (H): 

F  ■  KL  (CL  -  CS)  (2.3a) 

where  CS  ■  Pv/H  (2.3b) 

and   1  =  A  *  ET  (2.3c) 

KL   kL   HT£ 

in  which  <l  is  the  overall  mass  transfer  (or  liquid  film)  coefficient  based  on 
the  liquid  phase  (meters/day);  R  is  the  gas  constant  (atm.  m3/mol.K);  T  is 
the  temperature  (°K);  Pv  is  the  atmospheric  partial  pressure  (std.  atm.); 
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FIGURE  2.1  -  SCHEMATIC  REPRESENTATION  OF  TWO-FILM 
GAS  TRANSFER  MODEL 


and  Cs  1s  the  concentration  in  the  liquid  phase  in  equilibrium  with  Pv  (or 

saturation  concentration1.). 

■ 

The  resistances  to  mass  transfer  are  given  by  the  reciprocals  of  the  mass  transfer 
coefficients.  Therefore,  Eq.  2.3  can  be  written  as: 

rj  -  rL  +  rG  <2-4) 

where  rj  is  the  overall  resistance;  and  r|_  and  rG  are  the  resistances 
offered  by  the  liquid  ami  gas  films,  respectively.   Then,  the  fractions  of  the 
resistances  to  mass  transfer  in  the  liquid  and  gas  films  are  given  by 

a  =        i  <2-5a) 

rj   1  ■+  {rG/ri) 

IE  =     }  (2*5b) 

H"   1  +  (a/rG) 

where  rj_  =  _1  (2.5c) 

and  rG  =  RT/HkG  (2-5d) 

Equations  2.5a-b  can  be  utilized  to  determine  the  relative  magnitudes  of  the 
resistances  to  mass  transfer  in  the  liquid  and  gas  phases,  respectively.  The 
resistance  may  be  dominant  in  either  the  liquid  or  the  gas  phase,  or  in  both 
phases,  depending  on  th«  relative  magnitudes  of  k|_,  kG  and  H. 

From  Equation  2.5a,  the  smallest  *  of  resistance  in  liquid  film  is  seen  to  occur 
for  the  largest  kL  and  the  smallest  kg.  Using  the  reported  maximum  value 
of  kL  =  5.77  m/day  {from  the  data  on  benzene,  chloroform,  methylene  chloride 
and  toluene)  and  kG  =  480  m/day  (based  on  evaporation  of  water  in  a  canal),  and 
correcting  for  the  different  molecules  using  the  square  root  of  inverse  ratio  of 
molecular  weights,  the  %   resistance  in  liquid  film  for  ethylene  and  propane  were 
found  to  be  99.69*  and  99.91*,  respectively,  at  25°C  {Rathbun  and  Tai ,  1982). 
Based  on  this  analysis,  it  has  been  concluded  that  the  resistance  to  volatilization 
of  ethylene  and  propane  from  water  is  mostly  in  the  liquid  film. 

Rathbun  and  Tai  (1982)  have  also  evaluated  the  percentage  of  resistance  in  liquid 
film  as  a  function  of  the  Henry's  law  constant  (H)  by  utilizing  the  data  on  kL, 
kG  and  H  for  organic  compounds,  reported  in  the  literature  (see  Figure  2.2). 
Their  results  show  that,  for  average  conditions,  more  than  90*  of  the  resistance  is 





HENRY'S  LAW  CONSTANT  (atm-m3/mole) 


FIGURE  2.2  -  RESISTANCE  IN  LIQUID  FILM  AS  A  FUNCTION 
OF  HENRY'S  LAW  CONSTANT 
Source:  Rathbun  and  Tai (1982) 


1n  the  liquid  film  for  compounds  with  H  values  of  about  10"3  at  m.  m3/g 
mole.  Figure  2.2,  developed  by  Rathbun  and  Tai  (1982)  can  be  utilized  to  determine 
the  organic  compounds  which  have  volatilization  characteristics  similar  to  ethylene 
and  propane.  (Note:  The  Henry's  law  constants  are  discussed  further  in  Chapter 

V). 

2.2.2    Effects  of  Molecular  Properties: 

The  molecular  properties  of  chemical  compounds  which  affect  k|.     values     include 
molecular  weight,  molecular  diameter  and  diffusivities. 

The  relationship  of  molecular  weight  M  to  kL     (or     KL     for     compounds     with 
Henry's  law     constants     greater     than     10" 3     atm    m3/mol     under  average 
environmental   conditions)   can  be  expressed  as 

KLcc     M*0*5  (2-6) 

The  relationship  of  molecular  diameter  to  K|_  is  usually  expressed  as 

where  dm  refers  to  molecular  diameter.  Rathbun  and  Tai  (1982  and  references 
cited)  report  that  molecular  diameters  were  used  successfully  in  some  cases  to 
adjust  the  Kj_  for  a  tracer  gas  to  a  value  for  organics.  However,  information 
on  molecular  diameters  for  organic  compounds  of  interest  is  not  readily 
available;  and  hence,  molecular  diameters  are  often  estimated  through  the  use  of 
critical  volumes  and  the  longest  dimension  of  a  model  of  an  organic  compound.  The 
adjustment  of  KL  values  based  on  critical  volumes  is  described  in.  Chapter  IV. 

The  liquid  film  coefficient  has  also  been  related  to  the  molecular  diffusivity  of 
organic  compounds  in  water  by 

KL  -  ■£  (2-8) 

where  Dm  is  the  molecular  diffusion  coefficient  at  a  given  temperature  of  river 
water,   and  n  is  an  exponent.   Rathbun  and  Tai  (1981)  utilized  data   on 
volatilization  for  a  number  of  chlorinated  hydrocarbons  presented  by  Oil  ling  (1977) 
in  conjunction  with  Dm  values  calculated  from  an  equation  presented  by  Hayduk 
and  Laudie  (1974)  to  obtain  a  value  of  n  *  1.19;  however,  the  95*  confidence  limits 
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were  ±0.64.   Roberts  and  Dandliker  (1983)  found  n  =  0.66  for  six  chlorinated 
and  fluorinated  hydrocarbons.  Other  works  cited  by  Mackay  and  Yuen  (1983)  have 
found  values  of  n  to  be  0.5  and  0.67. 

Equations  2.6-2.8  can  provide  the  basis  for  calculating  the  values  of  K|_  for 
given  organic  compounds  from  the  known  values  for  tracers.  Since  M  and  d  are  basic 
properties  of  organic  compounds.  Equations  2.6  and  2.7  can  be  applied  for  adjusting 
the  Kl  values.  However,  the  molecular  diffusivity,  Dm  is  dependent  on  the 
properties  of  river  water  as  well  as  chemical  compound;  and  hence  Equation  2.8  may 
not  adequately  describe  the  effect  of  Dm  on  Kl.  Thus,  it  is  reasonable  to 
conclude  that  Equations  2.6  and  2.7  are  applicable  for  adjusting  the  Kl  values 
of  one  compound  to  another  whereas  the  effects  of  Dm  on  Kl  need  to  be 
considered  in  a  broader  context  in  the  analyses.  The  applicability  of  Equations 
2.6  and  2.7  is  dealt  with  further,  in  Chapter  IV. 

2.3   Review  of  Previous  Studies 

2.3.1  The  Rathbun-Tai  Relationships  for  Natural  Streams: 

Rathbun  and  Tai  (1982)  have  correlated  Kl  values  at  20° C  for  ethylene  and 
propane  with  hydraulic  data  from  eight  streams  and  rivers.  The  total  numbers  of 
<L  values  utilized  in  the  analyses  were  54  and  53  for  ethylene  and  propane, 
respectively.  Different  forms  of  predictive  equations  were  evaluated  using  the 
criterion  of  a  normalized  root-mean-square  (RMS)  error  defined  by: 


E  = 


100 

Ul"  V 

.   In  KL          f 

pred' 

A 

r 

n 

in  KL        T 
EXP 

(2.9) 


where  E  ■  RMS  error  (X) 

n  =  total  number  of  data  points 

EXP  and  PRED  =  subscripts  which  refer  to  experimental  and 
predicted  values,  respectively,  and 


fcn  Kl    ■  average  value  of  the  logarithms  of  the 
EXP         experimental  Kl  values 

The  KlexP  values  ranged  from  1.68  to  7.54  m/day  for  ethylene  and  1.13  to  6.94 
m/day  for  propane.  The  range  of  channel  hydraulic  parameters  were  as  follows: 


Discharge:    0.047  -  5.95  m5/s 
Depth:      0.101  -  0.555  m 
Velocity:    0.050  -  0.439  m/s 
Bed  Slope:    0.0538  -  0.631  m/km 

The  correlation  relationships  of  Rathbun  and  Tai  (1982)  are  as  follows: 

Ethylene:  KL  =  146  Z0030  (US)0**S*  (2.10) 

Propane:  KL  =  141  Z0'195  (US)0*515  (2.11) 


I 
I 
I 
I 
I 
I 

The  RMS  errors  (E)  associated  with  Equations  2.10  and  2.11  were  25.6%  and  34.8%,  14 
respectively.  Since  these  relationships  are  empirical  in  nature,  they  are  subject 
to  limitations  in  their  applicability  to  other  streams  and  rivers,  and  for  other  | 
organic  compounds. 

2.3.2  The  Mackay-Yuen  Laboratory  Studies: 

I 

Mackay  and  Yuen  (1983)  measured  volatilization  rates  of  11  organic  compounds  of 
varying  Henry's  law  constants  (5.57  x  10~*  to  5.18  x  10" 3  atm.  A 

m'/mol  at  20°C)  in  a  wind-wave  tank  in  the  laboratory.   The  experiments 
were  carried  out  at  six  different  windspeeds  in  the  range  5.96  -  13.2  m/s.  From  ■ 
these  studies,  the  following  equations  for  K[_  (m/day)  are  suggested:  ■ 

KL  =  0.0864  +  294.624  {UjJSzi),    (U^O.3)       (2.12)  | 

KL  =  0.0864  +  1244.16  (U*'2  /  JSqL)  ,  (li<0.03)     (2.13)  | 

where  U#  =  air  side  friction  velocity  (m/s),  and  ?| 


$CL  =  Schmidt  Number  =   absolute  viscosity 

diffusivity  x  density 


I 
I 
I 
I 
I 


The  value  of  U^  is  obtained  from 

U+  =  U10  (6.1  +  0.63UU)l/2  (2.14) 

where  U   1s  the  windspeed  at  a  height  of  10  meters  (m/s). 

1  o 

Equations  2.12  and  2.13  do  not  include  such  factors  as  depth  and  velocity  of 
water;  and  hence,  they  have  a  limited  scope  for  application  to  natural  streams  and 
rivers. 

2.3.3  Relationships  of  Southworth: 

Relationships  between  Kl,  properties  of  chemical  compounds,  windspeed  and  water 
depths  and  currents  have  also  been  presented  in  a  recent  publication  (Lyman  et  al , 
1982).  The  equations  for  K|_  (m/day)  developed  by  Southworth  (1979)  are: 

KL  =  5.64  U0-969  HT,    (Uw  <  1.9  m/s)  (2.15) 

KL  =  5.64    U°-969  /3T  .exp  [0.526  (Uw  -  1.9)], (1.9  <  Uwind  <  5  m/s)    (2.16) 
zo.s73  \JH 

where  M  =  molecular  weight  of  organic  compound  (g/mole),  and  Uw  =  windspeed 
(m/s).  While  these  relationships  include  depth,  velocity,  windspeed  and  molecular 
weight,  they  do  not  include  other  factors  such  as  slope,  viscosity  and  molecular 
diffusivity.  Therefore,  Equations  2.15  and  2.16  'are  also  subject  to   some 
limitations  in  estimating  liquid  film  coefficients  for  organic  compounds. 


11 


IH        DATA  COLLECTION  AND  ANALYSES 

3.1   Liquid  Film  Coefficients  from  Tracer  Studies 

A  modified  tracer  technique  developed  by  Rathbun  et  al  (1975)  was  utilized  by  the 
Ontario  Ministry  of  the  Environment  to  gather  data  in  shallow  streams  and  rivers. 
The  details  of  this  technique  have  been  presented  in  the  MOE  publication  Water 
Resources  Paper  No.  13  (Bacchus,  1981).  In  this  technique,  originally  developed 
for  direct  measurement  of  atmospheric  reaeration  -rates  in  streams,  ethylene  or 
propane  gas  and  a  conservative  material  (eg.  rhodamine  WT  dye)  are  utilized  as  the 
tracers.  The  technique  is  based  on  the  observation  that  the  ratio  between  the  rate 
of  desorption  of  a  tracer  gas  from  water  and  the  rate  of  absorption  of  oxygen  by 
the  same  water  remains  constant  regardless  of  the  mixing  conditions  (Tsivoglou, 
1967).  The  data  gathered  in  these  studies  are  also  found  to  be  useful  to  determine 
the  volatilization  rates  applicable  to  organic  compounds  (eg.  Rathbun  and  Tai,  1981 
and  1982). 

The  field  procedure  adopted  by  the  Ontario  Ministry  of  the  Environment  involves 
simultaneous  injection  of  ethylene  gas  and  rhodamine  WT  dye  solution  (as  a  line 
source  across  the  channel)  for  a  known  duration  of  time  at  the  head  of  a  study 
stretch  of  river  channel,  and  monitoring  for  dye  and  ethylene  gas  concentrations  as 
a  function  of  time  at  selected  downstream  stations  as  the  plume  passes  those 
stations.  These  data  are  utilized  to  calculate  the  volatilization  (or  desorption) 
coefficient  (Kv)  for  a  channel  reach  using  the  relationship: 

(Cg/Cd)US  \  (3-1) 


Kv,  =  1  £n 


(Cg/Cd)dS 


where  cg  and  c^  =  concentration  of  ethylene  gas  and  dye,  respectively; 
us  and  ds  -  subscripts  denoting  upstream  and  downstream  stations, 
respectively; 
t  =  travel  time  in  the  reach  (days) 
Kv  =  desorption  coefficient  of  ethylene  gas  (day"1) 

referred  to  base  e  at  the  average  river  water  temperature. 

The  liquid  film  coefficient  for  ethylene  gas  in  each  reach  is  then  determined  from 
the  equation  (Rathbun  and  Tai,  1981) 
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KL  =  Kv  Z  (3.2) 

where  Z  1s  the  average  depth  of  flow  in  the  channel  reach  (meters),  and  K|_  is 
the  liquid  film  coefficient  for  ethylene  gas  (meters/day).  The  K|_  values  will 
be  utilized  to  develop  relationships  of  K|_  to  bulk  flow  characteristics  of 
river  channels  and  properties  of  organic  compounds. 

3.2   Study  Areas  and  Hydraulic  Characteristics  of  River  Channels 

Data  from  the  following  streams  and  rivers,  gathered  by  the  Water  Resources  Branch, 
Ontario  Ministry  of  the  Environment,  were  analyzed  to  determine  the  liquid  film 
coefficients: 

•  Avon  River  below  the  Stratford  Water  Pollution  Control  Plant  (WPCP) 

•  Baden  Creek  below  Baden  WPCP 

•  Grand  River  near  Woolwich  Township,  Waterloo   County   (Close   to 
Winterborne) 

•  Nith  River  below  Wellesley  WPCP 

•  Speed  River  near  the  Township  of  Pus! inch  downstream  from  Guelph 

Figure  3.1  shows  the  layouts  of  the  various  study  segments  of  the  rivers  and 
streams.  Three  surveys  were  carried  out  in  the  Grand  and  Speed  rivers  and  one 
survey  was  conducted  in  each  of  the  other  streams. 

The  field  survey  dates  for  various  streams  and  the  salient  hydraulic 
characteristics  of  the  stream  channels  are  presented  in  Table  3.1.  Generally,  the 
stream  channels  are  very  shallow,  the  channel  bed  being  composed  mostly  of  rocks, 
pebbles  and  sand.  The  ranges  of  hydraulic  parameters  of  the  streams  are  as 
follows: 

Discharge  (Q):        0.0453  -  7.29  rn3/s 
Reach  length  (L):      127.5  -  2555.0  m 
Channel  width  (B) 


Average  depth  (Z) 
Mean  velocity  (U) 


3.11  -  38.6  m 

0.092  -  1.09  m 

0.028  -  0.27  m/s(2,419.2  -  23,328.0  m/day) 


Channel  bed  slope  (S):  0.069  -  0.0635  m/km  (0.00069  -  0.00635  m/m) 
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FIGURE  3.1   -  SCHEMATIC  DIAGRAMS  OF  STREAM  AND  RIVER  SEGMENTS 


TABLE  3.1 
SUMMARY  OF  CHANNEL  HYDRAULIC  DATA 


STREAM 

SURVEY 
DATE 

REACH 

L 
(m) 

Q 
(m3/s) 

B 

(m) 

7.31 

Z 
(ra) 

0.108 

U 
(m/d) 

12009.6 

S 

<m/m) 

AVON  R. 

9/04/80 

BC 

148.00 

0.110 

0.0034 

CD 

127.25 

0.110 

7.23 

0.092 

14256.0 

0.0034 

BADEN  CR. 

7/27/77 

BC 

274.00 

0.045 

3.11 

0.118 

10627.2 

0.0063 

CD 

1180.00 

0.045 

3.47 

0.236 

4752.0 

0-0063 

GRAND  R. 

6/02/78 

AC 

2555.00 

4.090 

29.47 

0.680 

17280.0 

0.0013 

GRAND  R. 

7/11/78 

AB 

115.00 

4.270 

29.69 

0.685 

18144.0 

0.0013 

BC 

2240.00 

4,270 

29.69 

0.685 

18144.0 

0.0013 

AC 

2555.00 

4.270 

29.69 

0.685 

18144.0 

0.0013 

GRAND  R. 

8/23/78 

BC 

2440.00 

7.290 

34.18 

0.790 

23328.0 

CO  013 

NITH  R. 

8/03/77 

A3 

303.00 

0.153 

17.91 

0.305 

2419.2 

0.0006 

BC 

260.00 

0.153 

13.80 

0.264 

3628.8 

O.QOOfi 

SPEED  R. 

6/07/78 

AB 

1480.00 

3.360 

30.68 

1.090 

8726.4 

0..0011 

BC 

1720.00 

3.360 

38.60 

0.440 

17280.0 

0.0011 

CD 

650.00 

3.360 

38.60 

0.440 

17280.0 

0 . 0O 11 

BD 

2380.00 

3.360 

38.60 

0.440 

17280.0 

0.0011 

SPEED  R. 

8/10/78 

AB 

1480.00 

2.120 

29.84 

0.890 

6998.4 

0.0011 

BC 

1720.00 

2.120 

38.60 

0.403 

12009.6 

0.0011 

SPEED  R. 

8/24/78 

AB 

1480.00 

2.430 

30.09 

0.943 

7430.4 

0.0011 

BC 

1720.00 

2.430 

38.60 

0.414 

13392.0 

0.0011 

CD 

650.00 

2.430 

38.60 

0.414 

13392.0 

0.0011 

BD 

2380.00 

2.430 

38.60 

0.414 

13392.0 

0.0011 
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An  evaluation  of  the  readily  available  channel  hydraulic  data  showed  that  the 
channel  widths,  depths  and  velocities  for  the  following  survey  conditions  were 
Identical  (see  Table  3.1): 

Grand  River,  June  11,  1978  -  Reaches  AB  and  BC 
Speed  River,  June  7,  1978  -  Reaches  BC  and  CD 
Speed  River,  August  24,  1978  -  Reaches  BC  and  CD 

Therefore,  in  these  three  cases,  the  two  successive  reaches  were  combined  into  one 
reach  (i.e.,  Grand  River  -  Reach  AC;  and  Speed  River  -  Reach  BD).  The  hydraulic 
characteristics  for  these  cases  are  given  in  Table  3.1  for  the  individual,  as  well 
as  the  combined  reaches. 

Data  on  channel  bed  slopes  were  obtained  from  topographic  maps.  However,  the 
slopes  for  individual  reaches  could  not  be  determined  from  these  maps.  In  the 
analyses,  the  drop  in  height  within  each  reach,  given  by  (reach  length  times  the 
slope),  will  be  utilized. 

3.3   Calculation  of  Liquid  Film  Coefficients  for  Ethylene  Gas 

The  data  analyses  were  carried  out  with  the  aid  of  the  computer  program  HGREDYET, 
documented  in  the  Water  Resources  Paper  No.  13  (Bacchus,  1981).  The  program  output 
includes  reaeration  rates  (Ka  per  day)  at  the  river  water  temperature,  T°C, 
and   at  20° C,  calculated  from  Equation  3.1.   These  were   converted   to 
volatilization  rates  (Kv),  using  the  relationship: 

Kv  =  0.89  Ka  (3.3) 

where  the  value  0.89  is  the  ratio  of  the  desorption  coefficient  for  ethylene  gas  to 
the  absorption  coefficient  for  oxygen  (Rathbun  et  al ,  1975).  The  liquid  film 
coefficients  (KjJ  at  the  average  instream  temperature  (T°C)  were  then 
calculated  from  Equation  3.2.  The  results  are  tabulated  in  Table  3.2.  The  values 
of  Kl  for  the  study  reaches  of  various  streams  and  rivers  are  seen  to  range 
from  0.97  to  8.07  meters/day.  The  river  water  temperatures  in  various  streams 
ranged  from  16°  to  25°C.  It  should  be  noted  that  the  data  for  the  Grand 
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River  -  June  11,  1978,  and  Speed  River  -  June  7,  and  August  24,  1978,  surveys 
Include  the  K[_  values  for  the  individual  reaches  as  welt  as  the  combined 
reaches  stated  in  the  previous  section. 

The  values  of  K(_  at  20CC  are  also  presented  in  Table  3-2\   These  are 
utilized  in  the  multiple  regression  analysis  described  in  the  next  chapter. 

3.4   Other  Data 

The  development  of  relationships  between  K|_  and  the   channel   hydraultc 
parameters,  presented  in  the  following  chapter,  requires  data  on  densities  and 
viscosities  of  stream  waters  as  well  as  molecular  diffusion  coefficient  of  ethylene 
gas  (Dm)  for  the  instream  temperature  conditions.  The  densities  (p)  and 
absolute  viscosities  (u)  were  obtained  from  the  CRC  Handbook  of  Chemistry  arrd 
Physics  (1980).  The  diffusion  coefficients  were  computed  from  the  following 
equation  (Hayduk  and  Laudie,  1974;  Lyman  et  al,  1982): 

13.26  x  10"5 
D  =  (3.4> 

m       1  •  Ik    u     0-  58  9 

M     VB 

where  Dm  is  the  molecular  diffusion  coefficient  of  ethylene  gas  (an2 /sec);;; 

u  is  the  absolute  viscosity  of  water  (centipoise  or  10" 2  g/cm  sec); 

and  Vg  is  the  molar  volume  of  ethylene  gas  computed  by  the  LeBas  method, 

outlined  in  Chapter  V,  (Vg  =  44.4  cm3/mol).   The  values  of  D^  (in 

cm2 /sec)   obtained  from  Equation  3.4,  have  been  converted'  to  the  units 

(m2/day).  The  values  of  density,  absolute  viscosity  and  molecular  diffusion 

coefficient  are  given  in  Table  3.2. 
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TABLE    3.2 
VISCOSITY, DENSITY, DIFFUSIVITY   AND   LIQUID   FILM  COEFFICIENTS 


STREAM 

SURVEY 
DATE 

REACH 

T 
<°C) 

U          P 
(Kg/m-day)  (Kg/m3) 

(m2/day) 

KL 

8T°C 
(m/day) 

K 

L 

@20°C 
(m/day) 

AVON  R. 

9/04/80 

BC 

20.7 

90.74 

998.0559 

0.000116 

2.66 

2.61 

CD 

20.7 

90.74 

998.0559 

0.000116 

2.39 

2.35 

BADEN  CR, 

7/27/77 

BC 

25.0 

76.90 

997.0448 

0.000140 

3.57 

3.17 

CD 

25.0 

76.90 

997.0448 

0.000140 

3.04 

2.70 

GRAND  R. 

6/02/78 

AC 

19.0 

88.73 

998.4052 

0.000119 

1.80 

1.84 

GRAND  R. 

7/11/78 

AB 

17.0 

93.40 

998.7748 

0.000112 

3.30 

8.92 

BC 

17.0 

93.40 

998.7748 

0.000112 

4.29 

4.60 

AC 

17.0 

93,40 

998.7748 

0.000112 

4.64 

4.98 

GRAND  R. 

8/23/78 

3C 

22.0 

82.50 

997.7704 

0.000129 

8.07 

7.69 

NITH  R. 

8/03/77 

AB 

25.0 

76.93 

997.0448 

0.000140 

2.63 

2.33 

BC 

25.0 

76.93 

997.0448 

0.000140 

2.26 

2.01 

SPEED  R. 

6/07/78 

AB 

16.0 

95.82 

998.9430 

0.000109 

0.97 

1.07 

• 

BC 

16.0 

95.82 

998.9430 

0.000109 

3.81 

4.19 

CD 

16.0 

95.82 

998.9430 

0.000109 

1.95 

2.14 

BD 

16.0 

95.82 

998.9430 

0.000109 

3.01 

3.31 

SPEED  R. 

8/10/78 

AB 

16.7 

94.12 

998.8252 

0.000111 

3.10 

3.35 

BC 

16.7 

94.12 

998.8252 

0.000111 

3.74 

4.05 

SPEED  R. 

8/24/78 

AB 

19.0 

88.73 

998.4052 

0.000119 

5.73 

5.87 

BC 

19.0 

88.73 

998.4052 

0.000119 

3.31 

3.39 

CD 

19.0 

88.73 

998.4052 

0.000119 

5.27 

5.39 

BD 

19.0 

88.73 

998.4052 

0.000119 

4.33 

4.44 
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IV         DEVELOPMENT  OF  PREDICTIVE  RELATIONSHIPS 

4.1  General  Approach 

In  order  to  develop  relationships  between  the  liquid  film  coefficient  (K[_)  and 
the  channel  hydraulic  parameters,  two  approaches  have  been  utilized.  The  first 
approach  is  a  multiple  regression  analysis  and  the  second  approach  is  a  dimensional 
analysis  using  the  Buckingham  n-theorem.  The  regression  analysis  is  mainly 
intended  as  an  attempt  to  identify  the  most  significant  channel  hydraulic 
parameters  affecting  K|_.  These  regression  analysis  results  will  likely  aid  in 
selecting  the  parameters  for  the  dimensional  analysis. 

4.2  Relationships  from  Multiple  Regression  Analysis 

For  multiple  regression  analyses,  the  functional  relationships  between  K[_  and 
the  channel  hydraulic  parameters  were  considered  to  be  of  the  general  forms: 

KL  =  a0  Bai  Za2  Ua3  (LS)a"  (4.1) 

KL  =  a0  Zaz  (US)as  (4.2) 

where       K|_  =  liquid  film  coefficient  of  ethylene  gas  at  20°C  (m/day) 
B  =  channel  width  (m) 
Z  ■  average  depth  (m) 
U  =  mean  velocity  (m/day) 
L  ■  length  of  reach  (m) 
S  =  channel  bed  slope  (m/m) 
(LS)  =  drop  in  height  in  reach  (m) 
a  to  a  ■  regression  coefficients 

OS 

It  should  be  noted  that  the  general  form  of  Equation  4.2  is  similar  to 
Equation  2.10  obtained  by  Rathbun  and  Tai  (1982). 

Equations  4.1  and  4.2  were  linearized  by  taking  logarithms  on  both  sides.  The 
computer  program  RLSEP,  available  in  the  IMSL  package,  was  utilized  to  obtain 
regression  models  using  a  forward  stepwise  algorithm.  The  results  of  the  multiple 
regression  analyses,  summarized  in  Table  4.1,  show  that  the  variation  explained  is 
in  the  range  6.8  -  22.4%,  whereas  the  residual  standard  deviations  expressed  as  a  % 
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TABLE  4.1 
SUMMARY  OF  MULTIPLE  REGRESSION  ANALYSES 


I 
I 


RUN 

NO. 

REGRESSION 
EQUATION 

REGRESSION 
COEFFICIENTS 

%  VARIATION 
EXPLAINED 

RESIDUAL  STD.   DEV. 
AS  A  %  OF  RESPONSE 
MEAN 

OVERALL! 
F-VALUE 

1 

Ki  =a0B 

a0s2.Q2 
ai=0.15 

6.8 

44.3 

0.94     ' 

2 

KL=a„Za2 

a  q=3.45 
a2=0.12 

3.7 

45.0 

0.50     1 

3 

KL=a0U33 

a o=0. 21 
a3s0.3 

14.4 

42.4 

2.20 

4 

KL-a0(LS)a- 

a o=2. 90 
a  .,=0.46 

13.2 

42.7 

2.00    1 

5 

KL=a0Ua3(LS)a* 

a  o=0.46 
a  3=0. 20 
a^O.10 

18.1 

43.2 

1.32    1 

6 

KL=a„Za2Ua3(LS)a* 

a  o=0.45 
a2=0.05 
a  3=0. 21 
a<,=0.08 

18.6 

45.0 

0.84    J 

7 

KL=a0Bai2a2Ua3(LS)aw 

a  o=0. 65 

ai=0.25 

a2=-0.21 

a3=0.06 

an=0.18 

22.4 

46.1 

0.72 

8 

KL=a0Za2(US)as 

a0s2.42 
a2=0.19 
as=0.15 

12.9 

44.5 

0.89 

of  response  mean  (i.e.,  %  variation  not  explained)  range  from  42.4  to  46. IX.  The 
latter  statistic  is  the  same  as  the  RMS  error  (Equation  2.9)  which  was  utilized  as 
the  criterion  for  selecting  the  regression  models  by  Rathbun  and  Tai  (1982).  The 
results  for  run  #7  (having  the  higher  %  variation  explained)  suggest  that  KL  is 
likely  a  function  of  B,  Z,  U  and  (LS).  However,  the  variation  not  explained  by  the 
regression  model  (Equation  4.1)  is  also  quite  high  (46.1%),  indicating  the  possible 
effect  of  other  factors  such  as  diffusivity  and  viscosity.  It  should  be  noted  that 
the  analyses  presented  in  this  section  are  referred  to  the  temperature  of  20°C. 

The  results  for  run  #8  show  the  fit  of  Equation  4.2  to  the  data.  These  results  are 
seen  to  differ  from  the  relationship  for  ethylene  (Equation  2.10)  obtained  by 
Rathbun  and  Tai  (1982). 

Table  4.1  also  gives  the  F-statistic  values  (indicative  of  significance  levels)  for 
various  regression  equations.  In  general,  the  values  are  lower  than  2.2  indicating 
that  they  are  significant  at  levels  of  about  15%  and  higher.  The  results  for  the 
runs  5-7  show  an  increase  in  the  %  variation  explained  as  the  number  of  dependent 
variables  is  increased;  however,  there  is  also  an  increase  in  the  variation  not 
explained,  accompanied  by  a  decrease  in  the  overall  F-values.  These  results 
suggest  a  somewhat  weak  statistical  correlation  of  K[_  to  the  channel  bulk  flow 
parameters. 

In  a  relative  sense,  the  parameters  U,  LS,  B  and  Z  are  important  in  that  order 
(based  on  the  overall  F-values  and  %  variation  explained).  Since  the  values  of 
slopes  are  taken  to  be  the  same  for  successive  reaches  for  some  of  the  streams 
(because  of  the  nature  of  the  data  available),  the  significance  of  the  term  (LS) 
may  not  be  discernible  from  this  analysis.  This  regression  analysis  suggests  that 
for  dimensional  analysis,  it  is  desirable  to  consider  all  the  four  parameters. 

4.3   Relationships  from  Dimensional  Analysis 

Dimensional  analysis  is  a  means  of  developing  relationships  between  dependent  and 
independent  variable  parameters  in  hydraulic  transport  studies.  It  involves  a 
systematic  organization  of  the  variable  parameters  into  the  smallest  number  of 
significant  dimensionless  parametric  groups.  Such  an  organization  is  facilitated 
by  the  fact  that  a  mathematical  relationship  must  be  dimensional ly  homogeneous;  in 
other  words,  each  term  in  an  equation,  when  reduced  to  the  basic  dimensions  of 
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length,  time  and  mass  (car  force),  must  contain  identical  powers  of  each  of  the 
respective  dimensions-  Thus,  if  the  variables  to  be  included  in  the  relationship 
are  known,  then  these  variables  should  occur  in  the  relationship  so  as  to  attain 
dimensional  homogenity.  The  principal  tool  of  dimensional  analysis  is  known  as  the 
Buckingham  n-theorem  (ftause,  1961). 

The  variable  parameters  which  can  be  considered  in  this  dimensional  analysis  are  as 
follows: 

liquid  film  coefficient  at  a  given  temeprature  of  river  water 

width,  depth,  velocity  and  drop  in  height  in  a  reach  of  stream  channel 

acceleration  dure  to  gravity 

density  of  water  at  the  given  temperature  of  river  water 

absolute  viscosity  of  water  at  the  water  temperature 

molecular  diffusivity  of  organic  compound  in  river  water 

wind  velocity 

The  molecular  weight  and  molecular  diameter  are  not  included  in  the  dimensional 
analysis  since  these  two  parameters  represent  the  properties  of  organic  compounds, 
as  stated  earlier  (see  Chapter  II).  Also,  data  on  wind  velocities  for  the  field 
survey  conditions  out1i?*ed  in  Chapter  III,  are  not  readily  available;  and  hence, 
wind  velocity  will  be  ignored  in  this  analysis.  It  should  be  noted  that  this 
analysis  is  specific  to  the  tracer  gas  under  consideration,  (i.e.,  ethylene  gas). 
However,  the  liquid  filra  coefficient,  density,  viscosity  and  diffusion  coefficient 
at  the  average  instreasn  temperature  will  be  considered  in  order  to  retain  the 
generality  of  the  relationship.  Thus,  the  general  relationship  can  be  written  as 

Kt  =  f(e,Z,U,LS,g,p,p,Dm)  (4.3) 

where  Kl  -  liquid  film  coefficient 

B  *  channel  width 

Z  =  average  depth 

U  -  mean  velocity 

LS  =  drop  in  height 

L  =  length  of  reach 

S  *  bed  slope 

g  =  acceleration  due  to  gravity 
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p  '*  density  of  water  at  T°C 
p  =  absolute  viscosity  of  water 
Dm  =  molecular  diffusivity  of  organic  compound 

Because  of  the  mathematical  equilibrium  between  the  independent  variable  (K|_) 
and  the  dependent  variables  appearing  on  the  right  hand  side  of  Equation  4.3,  the 
variables  can  be  grouped  in  another  functional  relationship  equal  to  zero: 

f  (KL,B,Z,U,lS,g,p,u,Dm)  =  0  (4.4) 

Then,  according  to  the  Buckingham  n-Theorem  (Rouse,  1961),  these  nine  variables 
can  be  grouped  into  (9-3)  =  6  dimensionless  n-terms  since  each  variable  can  be 
described  by  the  three  fundamental  dimensional  units  mass,  length  and  time  as 
follows: 

4>  =  (n  ,  n  ,  n  )  =  o  (4.5) 

Each  n-term  will  have  four  variables,  only  one  of  which  need  be  changed  from 
one  term  to  another. 

Numerous  combinations  of  the  parameters  into  the  six  dimensionless  n-terms  are 
possible.  For  this  analysis,  a  number  of  different  combinations  were  attempted. 
Correlations  of  the  data  to  these  relationships  resulted  in  varying  degrees  of 
variation  explained  (approximately  3-66*}.  Based  on  these  detailed  analyses,  the 
following  dimensionless  relationship  has  been  selected: 


W  sin^Hf  r($r(f  r 


(4.6) 


in  which  a0  is  a  numerical  constant  (dimensionless),  and  a^  -  as  are 
the  exponents  of  the  individual  dimensionless  groups.   The  ratio  of  absolute 
viscosity  to  density  (i.e.,  u/p)  is  equal  to  the  kinematic  viscosity 
(v).  Each  dimensionless  term  appearing  on  the  right  hand  side  is  related  to  a 
characteristic  property  affecting  K(_,  as  follows: 
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B  =  Aspect  ratio  of  channel 
Z 

LS  =  Ratio  of  drop  in  height  to  channel  depth 
1 

uZ  =  Rn  =  Reynolds  number  of  flow  indicative  of  viscous  effects 
v 

u2  =  f  =  Froude  number  denoting  gravitational  effects 

UZ  =  /  UZ  W  v  \=   Rn$cl  =  Product  of  Reynolds  (RN)  and 
0^   l7Mom)         Schmidt  (SC|_)  numbers 

Equation  4.6  is  linearized  by  taking  logarithms  on  both  sides.   The  values  of 
various  dimensionless  terms  were  calculated  by  using  the  data  summarized  in 
Tables  3.1  and  3.2  (see  Chapter  III).  The  IMSL  routine  RLSEP  was  utilized  to 
determine  the  regression  coefficients  a0  -  95.  From  this  procedure,  the 
following  relationship  is  obtained: 

IN 


The  various  statistical  parameters  associated  with  this  correlation  relationship 
are  as  given  below: 

Variation  explained  by  the  fit  =  66.6% 

RMS  error  (defined  by  Equation  2.9)  =  6.2% 

Overall  F-value  =  3.59 

Significance  level  =  4.6% 

These  statistical  characteristics  of  the  regression  fit  indicate  that  Equation  4.7 
satisfactorily  describes  the  functional  relationship  between  K[_  and  the  channel 
bulk  flow  parameters  and  chemical  characteristics  for  the  range  of  data  presented 
in  Tables  3.1  and  3.2. 

The  F-statistic  values  associated  with  individual  dimensionless  terms  were  in  the 
range  0.56  -  4.73.  The  F-values  for  the  last  three  terms  of  Equation  4.6  were 
significant  in  the  range  6  -  13%  whereas,  the  significance  levels  for  the  first  two 
terms  exceeded  20%.  Therefore,  a  regression  fit  was  carried  out  by  considering  the 
last  three  terms  of  Equation  4.6,  to  obtain  the  following  relationship: 
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Dm 
The  statistical  characteristics  of  the  regression  fit  are  as  follows: 

Variation  explained  =  62.14% 

RMS  error  (defined  by  Equation  2.9)  =  5.94% 

Overall  F-value  =  6.02 

Significance  level  =  l.K 

The  F-values  associated  with  the  individual  terms  were  in  the  range  2.44  -  6.18 
which  correspond  to  significance  levels  of  3-5%  approximately.  Thus,  Equation  4.8 
also  appears  to  describe  the  functional  relationship  satisfactorily.  This 
relationship  can  be  utilized  when  data  on  channel  slopes  are  not  readily  available. 

Equations  4.7  and  4.8  can  be  simplified  to  obtain  the  following  relationships: 

KL  =  2.70  g  °-3S  v_1-BS  Dm  1*89  B0'"  Z=0-2S  U0-26  (LS)0'16      (4.9) 


K,  =  7.16  g  °-23  v  -1"51  Dm1'59  Z0-15  U0-"6  (4.10) 

L 

where  v   =  u/p  is  the  kinematic  viscosity  of  water  at  the  instream  temperature 

According  to  Equation  4.9,  K|_  would  decrease  with  an  increase  in  depth  whereas 
Equation  4.10   indicates  K|_  increasing  with  depth.   Based  on  the  known 
relationships  of  reaeration  coefficient  to  channel  hydraulic  parameters,  the 
exchange  rates  at  the  air-water  interface  tend  to  decrease  with  increasing  depths 
(due  to  lower  turbulence  levels)  and  vice  versa.   Thus,  it  is  thought  that 
Equation  4.9  is  superior  to  Equation  4.10. 

Equations  4.7  to  4.10  are  dimensionally  homogeneous,  and  hence,  the  units  of  K[_ 
are  dependent  on  the  units  of  various  parameters  appearing  on  the  right  hand  side. 
For  example,  if  the  various  parameters  are  in  the  gram-meter-second  units,  then 
KL  is  in  meters/second.  Note  that  these  relationships  are  specifically  valid 
to  estimate  K^  for  ethylene  gas  at  a  given  temperature. 
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4.4   Evaluation  of  Predictive  Relationships 

In  this  section,  the  predictive  capabilities  of  Equations  4.9  and  4.10  will  be 
evaluated  by  a  comparison  with  the  observations  presented  in  Chapter  III  (for 
streams  and  rivers  in  Ontario)  and  the  data  presented  by  Rathbun  and  Tai  (1981). 
An  evaluation  of  the  applicability  of  Equation  2.10  of  Rathbun  and  Tai  (1982)  to 
the  provincial  streams  and  rivers  will  also  be  carried  out.  The  relationships  for 
K|_  developed  by  Mackay  and  Yuen  (1983)  and  Southworth  (1979),  presented  in 
Chapter  II,  could  not  be  evaluated  because  of  a  lack  of  appropriate  data. 

Equations  4.9  and  4.10  have  been  utilized  to  compute  the  values  of  K|_  for 
ethylene  gas  for  the  streams  and  rivers  listed  earlier  (see  Table  3.1).  The 
computed  values  of  K|_  are  presented  in  Table  4.2  along  with  the  observed  K|_ 
values  (range  of  observed  K(_:  0.97  -  8.07  m/day).  The  values  obtained  from 
Equation  4.9  are  in  the  range  1.87  -  5.57  m/day,  whereas  Equation  4.10  gives  values 
in  the  range  of  2.51  -  6.31  m/day.  Figures  4.1  and  4.2  show  plots  of  the  K^ 
values  computed  from  Equations  4.9  and  4.10  versus  the  observed  K|_  values.  The 
average  K[_  values  are  3.46  m/day  for  the  observations,  and  3.03  and  3.28  m/day 
from  Equations  4.9  and  4.10,  respectively.  The  ratios  (predicted  Ki_/observed 
K(J  vary  from  0.44  to  2.11  with  Equation  4.9,  and  0.51  to  2.64  with  Equation 
4.10;  and  the  average  ratios  are  1.01  and  1.12  for  Equations  4.9  and  4.10, 
respectively.  These  results  indicate  that  the  predictions  of  Equations  4.9  and 
4.10  compare  favourably  with  the  observations.  On  an  overall  basis,  Equation  4.9 
is  slightly  better  than  Equation  4.10. 

The  values  of  K[_  for  ethylene  gas,  computed  from  Equation  2.10  of  Rathbun  and 
Tai  (1982),  are  presented  in  Table  4.2.  The  computed  values  (K|_rj)  are  in  the 
range  0.74  -  3.44  m/day.  Figure  4.3  shows  a  plot  of  the  computed  versus  observed 
K(_  values.   Generally,  the  computed  values  are  seen  to  be  lower  than  the 
observations;  the  averages  of  computed  and  observed  values  are  2.1  and  3.46  m/day, 
respectively.  The  predicted  values  range  from  0.34  to  2.27  times  the  observed 
values,  the  average  ratio  being  0.75  (see  Table  4.2).  These  results  show  that 
Equation  2.10  developed  by  Rathbun  and  Tai  (1982)  generally  underestimates  the 
K|_  values  for  the  provincial  streams  and  rivers. 
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TABLE  4.2 
COMPARISON  OF  OBSERVED  AND  COMPUTED  K^  FOR  ONTARIO  STREAMS 


OBSERVED 

\ 

(m/day) 

EQUATION  4.9 
(m/day)   K 

EQUATION  4.10 
(m/day)   K 

EQUATION  2.10 
(m/day)   K 

2.66 

2.10 

0.79 

2.51 

0.95 

2.01 

0.75 

2.39 

2.22 

0.93 

2.66 

1.11 

2.06 

0.86 

3.57 

3.34 

0.94 

4.17 

1.17 

2.61 

0.73 

3.04 

3.01 

0.99 

3.19 

1.05 

2.26 

0.74 

1.80 

3.80 

2.11 

4.22 

2.35 

2.85 

1.58 

4.64 

3.14 

0.68 

3.65 

0.79 

2.92 

0.63 

8.07 

5.57 

0.69 

6.31 

0.78 

3.44 

0.43 

2.63 

2.57 

0.98 

2.43 

0.92 

0.64 

0.24 

2.26 

2.60 

1.15 

2.86 

1.27 

0.74 

0.33 

0.97 

1.87 

1.93 

2.57 

2.64 

2.20 

2.27 

3.01 

3.32 

1.10 

3.07 

1.02 

2.24 

0.74 

3.10 

1.97 

0.64 

2.39 

0.77 

1.86 

0.60 

3.74 

3.15 

0.84 

2.72 

0.73 

1.84 

0.49 

5.73 

2.51 

0.44 

3.01 

0.52 

1.94 

0.34 

4.33 

4.28 

0.99 

3.48 

0.80 

1.95 

0.45 

AVERAGE:   3.46 

3.03 

1.01 

3.28 

1.12 

2.10 

0.75 

KlC  :  COMPUTED  VALUE 
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FIGURE  4.1 
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OBSERVED  KL  (m/day) 

COMPARISON  OF  KL  FROM  EQUATION  4.  9 
WITH  THE  OBSERVED  VALUES 
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FIGURE  4.2 


COMPARISON  OF  KL  FROM  EQUATION  4.10 
WITH  THE  OBSERVED  VALUES 
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FIGURE  4.3 


3    4    5    6    7    8    9 

OBSERVED  KL  (m/day) 

COMPARISON  OF  KL  FROM  EQUATION  2.10  OF 
RATHBUN  AND  TAI  (1981)  WITH  THE  OBSERVED 
VALUES 


Rathbun  and  Tai  (1981)  have  presented  values  of  liquid  film  coefficient  at  20°C 
(Klq)  applicable  to  the  organic  compounds  benzene,  chloroform,  methylene 
chloride  and  toluene  for  several  streams  and  rivers.  The  values  of  K|_q  for 
these  four  compounds  were  found  to  be  approximately  the  same,  based  on  laboratory 
studies.  These  Klq  values  are  related  to  the  Kl  values  for  ethylene  by  the 
relationship  Kj_q  -   0.753  Kl.  The  river  hydraulic  data  given  by  Rathbun  and 
Tai  (1981)  include  water  depth,  velocity  and  discharge;  but  the  bed  slopes  and 
reach  lengths  are  not  reported.  These  data  were  thus  utilized  to  compute  the 
values  of  liquid  film  coefficients  for  ethylene  from  Equation  4.10.  Table  4.3 
shows  the  hydraulic  data  for  the  streams  as  well  as  the  K|_  values  obtained  from 
KLq  ■  0-753  KL  (Rathbun  and  Tai,  1981).  Herein,  these  calculated  KL 
values  will  be  referred  to  as  "observed"  values.   The  range  of  the  various 
parameters  are  as  follows: 

Discharge  0.027  -  93.50  m3/s 

Water  depth  0.053  -  2.330  m 

Velocity  0.060  -  0.747  m/s 

K|_  0.46  -  7.66  m/day 

The  predictions  of  Equation  4.10,  given  in  Table  4.3,  are  seen  to  range  from  1.87 
to  10.10  m/day.  The  average  of  the  "observed"  and  computed  K|_  values  are  2.00 
and  4.69  m/day,  respectively;  and  the  ratios  between  the  "computed"  and  observed 
values  range  from  0.75  to  7.79,  the  average  ratio  being  3.53.  These  results  show 
that  Equation  4.10  overestimates  K|_  values  for  the  streams  and  rivers  for  which 
data  are  reported  by  Rathbun  and  Tai  (1981).  It  is  recalled  that  Equation  2-10 
developed  by  Rathbun  and  Tai  (1982)  underestimates  K[_  values  for  provincial 
streams  and  rivers.  These  discrepancies  are  likely  to  be  due  to  the  possible 
effects  of  such  factors  as  wind  velocity,  bed  slope,  and  the  presence  of  other 
chemical  compounds  in  the  streamwaters.  It  is  also  possible  that  the  indirect 
method  of  estimating  the  "observed"  K|_  values  for  ethylene  has  contributed,  to 
some  extent,  for  these  discrepancies. 
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TABLE  4.3 

COMPARISON  OF  PREDICTIONS  OF  EQUATION  4.10 
WITH  DATA  OF  RATHBUN  AND  TAI  (1981) 


z 

U 

"lo 

K 
LC 

^C 

(m) 

(ra/s) 

(m/day) 

(m/day) 

\o 

0.472 

0.160 

2.31 

3.91 

1.69 

0.285 

0.305 

2.51 

4.88 

1.94 

0.289 

0.439 

7.66 

5.78 

0.75 

0.387 

0.317 

2.47 

5.20 

2.11 

0.053 

0.160 

2.12 

2.82 

1.33 

0.090 

0.150 

1.46 

2.96 

2.03 

0.082 

0.150 

1.33 

2.92 

2.20 

0.070 

0.060 

0.46 

1.87 

4.03 

0.260 

0.160 

0.88 

3.58 

4.08 

0.270 

0.085 

0.48 

2.69 

5.63 

0.180 

0.210 

0.92 

3.84 

4.19 

0.280 

0.180 

1.73 

3.82 

2.21 

0.482 

0.160 

1.04 

3.92 

3.77 

0.290 

0.270 

1.39 

4.63 

3.32 

0.518 

0.408 

0.78 

6.10 

7.79 

0.579 

0.552 

5.79 

7.13 

1.23 

0.270 

0.088 

0.68 

2.73 

4.04 

0.948 

0.144 

1.25 

4.14 

3.32 

0.701 

0.171 

4.78 

4.28 

0.90 

1.240 

0.546 

0.81 

7.95 

9.82 

1.270 

0.539 

.  1.46 

7.94 

5.43 

2.330 

0.747 

1.73 

10.10 

5.85 

AVERAGE: 

2.00 

4.69 

3.53 

KQ  :" OBSERVED "  ;     &r  : COMPUTED 
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4.5   Sensitivity  Analysis 
4.5.1  General  Procedure: 

A  knowledge  of  the  sensitivity  of  various  dependent  parameters  on  the  independent 
variables  appearing  in  Equations  4.7  -  4.10  will  aid  in  identifying  the  relative 
importance  of  the  dependent  parameters  in  estimating  K^  Herein,  true 
sensitivity  analyses  will  be  carried  out  by  the  method  of  determining  the  relative 
error  or  relative  change  (described  in  standard  text  books  on  calculus). 
Equations  4.7  -  4.10  are  of  the  general  form 

Y  =  a0  X^1  X*2  Xjjr.  C4U1), 

where  Y  is  the  independent  variable; 

Xi  (1=1, n),  are  the  dependent  variables; 
n  is  the  number  of  dependent  variables; 
a0  is  a  constant;  and 

&\   is  the  exponent  of  X-j  (i=l,n) 

In  order  to  determine  the  relative  change  in  Y  due  to  a  given  relative  chance  in 
%l,   Equation  4.11  is  differentiated  with  respect  to  Xj  to  obtain 

^•a0alXlX2!  Xn"  S4-I2) 

From  Equations  4.11  and  4.12,  the  following  relationship  is  obtained; 

-£*a,dXl  (4.13) 

In  Equation  4.13,  (dXi/Xi)  and  (dY/Y)  denote  relative  changes,  fn  X^  and 
Y,  respectively.  Thus,  for  a  given  relative  change  in  X^,  the  relative  change 
in  Y  can  be  obtained  from  Equation  4.13.  In  other  words,  the  sensitivity  of   iJfj 
can  be  determined  from  this  relationship.  Similar  analyses  can  be  carried  out  for 
the  other  variables  X,  (i=2,n)  appearing  in  Equation  4.11. 
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4.5.2  Sensitivity  of  Various  Parameters: 

The  general  procedure  described  in  the  previous  section  has  been  applied  to  examine 
the  sensitivities  of  various  parameters  of  Equations  4.7  -  4.10.  The  results  of 
sensitivity  analyses  are  presented  in  Table  4.4.  In  these  analyses,  the  relative 
change  in  the  independent  variable  due  to  a  relative  change  of  20%  in  the  dependent 
variable  under  consideration,  has  been  determined.  The  results  shown  in  Table  4.4 
are  arranged  in  a  decreasing  order  of  relative  importance  of  various  dependent 
parameters.  The  most  sensitive  parameters  are  the  molecular  diffusivity  of  organic 
compound  and  kinematic  viscosity  of  water. 

4.6   Relation  Between  Liquid  Film  Coefficients  of 
Tracer  Gases  and  Organic  Compounds 

4.6.1  Review  of  Relationships: 

Based  on  studies  carried  out  by  Rathbun,  et  al  (1978)  and  Rainwater  and  Hoi  ley 
(1984),  the  liquid  film  coefficients  of  any  two  volatile  compounds  are  related  by 

R  -  !^_  (4.14) 

klorg 

where  KLG  and  KL0RG  are  the  liquid  film  coefficients  of  a  tracer  gas  and  organic 
compound,  respectively;  and  R  is  the  ratio  between  the  two  coefficients.  Rathbun, 
et  al  (1978)  and  Rainwater  and  Hoi  ley  (1984)  have  presented  values  of  R  between 
various  tracer  gases,  determined  from  laboratory  experimental  studies.  Also,  the 
value  of  R  between  ethylene  gas  and  a  group  of  four  organic  compounds,  has  been 
given  by  Rathbun  and  Tai  (1981)  as  stated  earlier  (see  Section  4.4). 

The  values  of  R  can  also  be  calculated  by  making  use  of  molecular  diameter, 
critical  volume  or  molecular  weight  of  organic  compounds  (see  Chapter  II).  In  a 
recent  publication,  Rainwater  and  Holley  (1984)  have  presented  theoretical  aspects 
related  to  the  use  of  molecular  diameter  and  critical  volume  in  determining  the 
values  of  R.  The  relationship  obtained  by  Rainwater  and  Holley  (1984)  can  be 
written  in  the  form 
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TABLE  4.4 
RESULTS  OF  SENSITIVITY  ANALYSES 


Dependent  Variable 

Independent  Var 

iable 

Relationship 

Parameter 

Relative  Change 

Relative  Change 

Parameter 

Equation  4.7 

UZ/Dm 

20.0 

-37.8 

UZ/v 

20.0 

37.0 

B/Z 

20.0 

8.0 

KL/U 

uVgZ 

20.0 

7.0 

LS/Z 

20.0 

3.2 

Equation  4.8 

uz/om 

20.0 

-31.8 

UZ/v 

20.0 

30.2 

KL/U 

u2/gz 

20.0 

-  4.6 

Equation  4.9 

Dm 

20.0 

37.8 

V 

20.0 

-37.0 

B 

20.0 

8.0 

KL 

U 

20.0 

5.2 

z 

20.0 

-  5.0 

LS 

20.0 

3.2 

Equation  4.10 

Dm 

20.0 

31.8 

V 

20.0 

-30.2 

«L 

u 

20.0 

9.2 

Z 

20.0 

3.0 
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where  Rv  is  the  ratio  based  on  critical  volumes;  dG  and  dgRG  denoted, 
1n  order,  the  molecular  diameters  of  tracer  gas  and  organic  compound;  and  Vqg 
and  Vcorg  are  the  critical  volumes  of  the  tracer  gas  and  organic  compound, 
respectively.  This  relationship  is  based  on  the  assumption  that  the  molecules  are 
spherical,  and  hence,  the  critical  volume  is  proportional  to  the  cube  of  the 
diameter.  The  exponent  0.942  appearing  in  Equation  4.15  is  very  close  to  unity, 
and  hence,  the  relationship  between  R  and  the  ratio  of  molecular  diameters  is  a 
reasonable  agreement  with  Equation  2.7  (i.e.,  K[_   dm  ). 

Based  on  the  relationship  of  Kj_  to  molecular  weight  presented  in  Chapter  II 
(Equation  2.6),  the  following  expression  is  obtained. 

R   _   KLG   /MqRG  \0'5  (4-16) 

KLORG   \  Mq   / 


where  Rh  is  the  ratio  based  on  the  molecular  weights,  and  Mq  and  Mqrg 
denote,  in  order,  the  molecular  weights  of  tracer  gas  and  organic  compound.   A 
rigorous  evaluation  of  Equation  4.16  (similar  to  that  for  Equation  4.15)  could  not 
be  found  in  the  literature,  except  for  the  studies  of  Southworth  (1979)  described 
in  Chapter  II. 

This  review  indicates  that  either  one  or  both  of  the  above  relationships  (i.e., 

Equations  4.15  and  4.16)  can  be  the  basis  for  adjusting  the  K[_  values  of  tracer 

gases  to  organic  compounds  of  interest.  A  comparitive  study  of  observed  and 

computed  values  of  P.  will  aid  in  deciding  the  validity  of  Equations  4.15  and  4.16. 

4.6.2  Evaluation  of  Relationships: 

Laboratory  experimental  data  presented  by  Rathbun,  et  al  (1978),  Rathbun  and  Tai 
(1981),  Roberts  and  Dandliker  (1983),  Kaczmar,  et  al  (1984)  and  Mackay  and  Yuen 
(1983)  will  be  utilized  to  evaluate  the  validity  of  Equations  4.15  and  4.16. 
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Rathbun  and  Tai  (1981)  present  K(_  values  for  oxygen  and  four  organic  compounds 
whereas  Rathbun  et  al  (1978)  report  desorption  -  absorption  rates  (per  day) 
associated  with  ethylene,  propane  and  oxygen.  The  studies  of  Roberts  and  Dandliker 
(1983),  and  Kaczmar,  et  al  (1984)  include  absorption  rates  of  oxygen  (per  day)  and 
desorption  rates  (termed  mass  transfer  rate  constants  expressed  in  the  units  per 
second)  for  several  volatile  compounds.  Thus,  exchange  rates  for  oxygen  are 
measured  in  all  of  these  studies.  The  analyses  presented  below  involve  ratios  of 
exchange  rates  between  oxygen  and  an  organic  compound.  These  ratios  have  been 
determined  by  using  either  the  K|_  values  or  the  absorption  and  desorption 
rates.  For  the  sake  of  convenience,  the  exchange  rates  (i.e.,  liquid  film 
coefficients,  and  absorption  and  desorption  rates)  will  be  denoted  by  K'. 

The  following  relationships  permit  computation  of  the   ratios   of   exchange 
coefficients: 

Rf   ■  (K1  for  oxygen/K'  for  ethylene)  (4.17) 

R02  ■  (K'  f°r  oxygen/K'  for  organic  compound)  (4.18) 

RORG  "   (Roz/RE)  =  (K'  for  ethylene/*'  for  organic  compound)   (4.19) 

The  ratios  obtained  from  Equations  4.17  -4.19  will  be  termed  the  observed  ratio 
values  since  they  are  obtained  from  measured  K'  values.  Table  4.5  shows  the 
observed  Rq2  and  Rqrg  values  for  various  volatile  compounds  including 
oxygen.  The  Henry's  law  constants  and  molecular  weights  of  the  compounds  are  also 
given  in  this  table. 

Mackay  and  Yuen  (1983)  obtained  liquid  film  coefficients  (KL)  at  six  wind 
speeds  (range  5.96  -  13.2  m/s)  for  eleven  organic  compounds.  The  compounds  tested 
as  well  as  their  Henry's  law  constants,  molecular  weights  and  critical  volumes  are 
given  in  Table  4.5.   Since  oxygen  has  not  been  included  in  this  study,  the  values 
of  Rqrg,  defined  by  Equation  4.19,  were  calculated  indirectly.   First,  the 
ratios  (Rbz)  were  computed  by  considering  benzene  as  the  reference  compound. 
Then,  these  ratios  were  converted  to  the  Rqrg  values  through  the  use  of  the 
observed  ratio  (K"  for  ethylene/K'  for  benzene)  =  1.29.  The  relationships  utilized 
are: 
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TABLE  4.5 
RATIOS  OF  EXCHANGE  COEFFICIENTS  FOR  VOLATILE  COMPOUNDS 


Source 

or 

(lata 


Compound 


Henry's  Law  Constant 
(atm-m'/mol  at  20°C 

and  1  atm) 


Molecular   Critical 
Weight    Volume 
(gram/no  I)   (cm'/iwl) 


Observed  Ratio  *     Computed  Ratio 

Equation  4.15         Equation  4.16 
"02       R0RG       RV  02       RV0RG       RMQ2      BH0RG 


Rathbun  «t  al 


Rathbun  I   Tat 
(1981)° 


Roberts  and 

Oandllker 

(1983)c 


Kaczmar  et  al 
£    (1904)° 

Hackay  and  Yuen 
(1983)e 


Henry's  Law  Constant 
>  10*1  atni-m'/mol 

Oxygen 
Propane 

Benzene 
Toluene 

Methylene  chloride 
Chloroform 

Chloroform 

Carbon  tetrachloride 

Otchlorodlfluoromethane 

Tetrachloroethene 

1,1,1-trlchloroethane 

Trichloroethene 

Chlo  oform 
Bromodlchlorome thane 

Benzene 

Toluene 

1,2-dlchloropropane 

Chlorobenzene 

Carbon  tetrachloride 


0.73 

1.1  X  10  ' 

4.4  x  10"J 

5.2  -  6.6  x  10"' 
3.0  x  10"1 

3.4  x  10"' 

3.4  x   10"' 

2.3  -  2.5  x  10"' 

1.5 

8.3  x  10"* 

3.6  -  18  x  10"' 

1.0  x  10"' 

3.4  x  10"' 

1.3  x  10"' 

4.4  x   10'* 

5.2  -  6.6  x   10*' 

2.1  x   10   ' 
2.6  x  10"1 

2.3  -  2.5  x  10"* 


32.0 

73.4 

?fl;[ 

!20,Q 

4*J 

ffli.d 

78.1 

262.0 

92.2 

316.0 

84.9 

193.0 

119.4 

238.0 

119.4 

238.0 

153.8 

277.0 

120.9 

215.0 

165.8 

308.0 

133.4 

283.0 

131.4 

268.0 

119.4 

238.0 

163.8 

259.0 

78.1 

262.0 

92.2 

316.0 

113.0 

299.0 

112.6 

310.0 

153.8 

277.0 

1.00 

1:8 

1.49 
1.56 
1.45 
1.50 

1.79 
1.62 
1.54 
1.66 
1.67 
1.63 


1.59 
1.82 


I  :«8 

l.H 

1.29 
1.36 
1.26 
1.30 

1.56 

1.41 
1.34 
1.44 
1.45 
1.42 

1.38 
1.58 

1.29 
1.25 
1.33 
1.22 
1.24 


1.00 

rfi 

i.m 

1.49 
1.58 
1.35 
1.45 

1.45 
1.52 
1.40 
1.57 
1.53 
1.50 

1.45 
1.49 


Lis 

1.27 
1.35 

1.15 
1.24 

1.24 

1.30 
1.20 
1.34 
1.31 
1.28 

1.24 
1.27 

1.49 
1.58 
1.33 
1.35 
1.30 


1.00 

0.943 

1.17 

1:8 

1.56 

1.65 

1.70 

1.80 

1.63 

1.73 

1.93 

2.05 

1.93 

2.05 

2.19 

2.32 

1.94 

2.C6 

2.28  ■ 

2.42 

2.04 

2.16 

2.02 

2.14 

1.93 

2.05 

2.26 

2.40 

1.67 

1.80 

2.00 

2.00 

1.75 

Henry^s  law  Constant 
<  10'*  atm-m'/mol 


Kaczciar  et  al 
(1984)d 

Hack  ay  and  Yuen 
(1983)f 


Chlorodlbromomethane 

Bromoform 

1,2  dibroroo«ethane 

2-pentanone 

2-heptanone 

1-pentanol 

2-methyl-l-propanol 

l-butanol 


8.1  x 
5.8  x 


6.3  x  10"* 

3.2  x  10"' 

9.0  x  10"s 

1.0  x  10"' 

1.0  x  10" * 

5.6  x  10"' 


Average  of  the  number  of  tests  stated  below: 

?  -  81-  tests  for  ethylene  and  34  tests  for  propane 
"  -  TO  tests  for  each  of  the  four  compounds 
c  -  16  tests  for  each  of  the  six  compounds 


208.3 

280.0 

252.8 

301.0 

173.9 

290.0 

86.1 

320.0 

114.2 

430.0 

88.2 

333.0 

74.1 

274.0 

74.1 

278.0 

d 

e 

f 

2.38 

2.07 

3.45 

3.00 

1.56 

3.26 

2.77 

12.40 

14.40 

21.52 

1.52 

1.30 

1.56 

1.33 

1.32 

1.36 

1.49 

1.38 

1.30 

1.30 

2.55 

2.70 

2.81 

2.98 

2.49 

1.75 

2.02 

1.77 

1.62 

1.62 

-  13  tests  for  chloroform  and  2  tests  for  each  of  the  other  three  compounds 

-  6  tests  for  benzene  and  toluene,  4  tests  for  the  next  three  compounds 

-  4  tests  for  the  first  three  and  3  tests  for  the  last  three  compounds 


RBZ  =  (KL  for  benzene/KL  for  organic  compound)      (4.20) 
Rqrg  =  1.29  RBZ  (4.21) 

The  values  of  Rg7_  utilized  herein  are  the  average  values  of  the  ratios  obtained 
for  each  of  the  windspeeds.  Benzene  and  toluene  were  included  in  all  the  six 
windspeed  experiments,  whereas  the  remaining  compounds  were  monitored  at  three  or 
four  windspeeds.  The  observed  ratios  (Rqrg).  obtained  from  Equations  4.20  and 
4.21,  are  presented  in  Table  4.5. 

Four  of  the  compounds  listed  in  Table  4.5,  namely  benzene,  toluene,  chloroform  and 
carbon  tetrachloride,  were  utilized  in  two  or  more  studies.  The  observed  values  of 
RQRG  f°r  each  compound,  obtained  from  these  studies,  are  seen  to  differ 
somewhat.  This  indicates  the  possible  effects  of  testing  conditions  on  the  observed 
R0RG  values. 

The  critical  volumes  (Vr,)  for  various  volatile  compounds  were  determined  as 
outlined  in  Chapter  V  (see  Section  5.5).   The  compounds  are  listed  in  Table  4.5 
along  with  the  corresponding  values  of  (VrJ  and  molecular  weights  (M).   These 
data  have  been  utilized  to  calculate  the  values  of  R  and  R^  according  to 
Equations  4.15  and  4.16,  respectively.  The  computed  values  of  Ry  and  R^ 
for  various  compounds  are  given  in  Table  4.5. 

The  range  of  various  parameters  for  the  volatile  compounds  presented  in  Table  4.5 
are  as  follows: 

Molecular  weight:  28.1  -  252.8  grams /mo 1 

Critical  volume:  73.4  -  316.0  cm3/mol 

Observed  ratio  (R02)  between  oxygen  and:    1.00  -  3.45 
organic  compound 

Ratio  (Rqrg)  between  ethylene  and  organic:   1.00  -  21.52 
compound,  based  on  observations 

Ratio  (R\/02)»  computed  from  critical:       l.CO  -  1.58 
volumes  using  Equation  4.15 

Ratio  (R\/0RG)»  based  on  critical  volumes:    1.00  -  1.35 

Ratio  (Rm02)»  computed  from  molecular:      0.943  -  2.81 
weights  using  Equation  4.16 


Ratio  (Rmorg)»  based  on  molecular  weights:   1.00  -  2.98 

A  comparison  of  the  computed  Ry  and  Rm  values  with  the  corresponding 
observed  ratios  for  various  compounds  indicates  that: 

(1)  the  computed  values  of  Ry,  obtained  from  Equation  4.15  using  the  critical 
volumes,  are  in  reasonable  agreement  with  the  corresponding  observed  ratios 
with  the  exception  of  the  compounds  for  which  the  Henry's  law  constants  are 
smaller  than  10" 3  atm-m3/mol. 

(2)  the  computed  values  of  R^,  obtained  from  Equation  4.16  using  the  molecular 
weights,  differ  from  the  corresponding  observed  ratios  for  about  50%     of  the 
compounds  listed  in  Table  4.5;   however,  for  bromoform,  the  computed  Rm 
value  is  in  agreement  with  the  observed  ratio. 

(3)  the  observed  P.  values  for  chloroform  are  in  the  range  1.30  -  1.56,  indicating 
possible  effects  of  variations  in  the  testing  conditions. 

This  analysis  suggests  that  Equation  4.15  can  be  utilized  to  compute  the  values  of 
R  required  for  adjusting  Kj_  values  of  highly  volatile  compounds  (Henry's  law 
constant  >10"3  atm-m3/mol),  although  some  limitations  are  evident.  It 
1s  recommended  to  evaluate  the  validity  of  Equations  4.15  and  4.16  for  a  broader 
range  of  organic  compounds  as  data  become  available,  as  well  to  undertake 
Investigations  for  developing  similar  relationships  for  compounds  of  moderate  and 
low  volatility. 

4.7   Computation  of  Liquid  Film 

Coefficient  of  an  Organic  Compound 

The  various  steps  involved  in  computing  the  liquid  film  coefficient  for  a  given 
organic  compound  are  outlined  below.  The  computational  procedure  is  then 
illustrated  by  an  example. 
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Computational  Procedure: 

1.  For  the  given  organic  compound,  obtain  the  Henry's  law  constant  (see 
Chapter  V).  This  computational  procedure  is  generally  valid  for  compounds 
having  Henry's  constants  greater  than  10   atm  -  nr/mol. 

2.  Obtain  hydraulic  data  for  the  desired  reaches  (or  segments)  of  a  given  river 
channel.  These  data  include:  discharge,  reach  length,  width,  depth,  velocity 
and  bed  slope.  (Note:  If  the  latter  is  unknown,  Equation  4.10  could  be 
utilized.) 

3.  For  the  desired  river  water  temperature,  obtain  the  values  of  density  and 
absolute  viscosity,  which  are  readily  available  in  standard  fluid  mechanics 
text  books  or  handbooks  on  physical  and  chemical  properties  (eg.  Weast  and 
Astle,  1979). 

4.  Determine  the  molar  volume  (Vg)  for  ethylene  gas  by  the  LeBas  method  as 
outlined  in  Section  5.5. 

5.  Calculate  Dm  for  ethylene  gas  by  using  the  Hayduk-Laudie  relationship 
(Equation  3.4). 

6.  Compute  the  liquid  film  coefficient  (K|_E)   for   ethylene   gas   from 
Equation  4.9  or  4.10  (depending  on  whether  the  bed  slope  is  known  or  unknown). 

7.  Calculate  the  critical  volume  (cm3/mol)  for  the  organic  compound  as 
outlined  in  Section  5.5  (see  Chapter  V). 

8.  Determine  the  ratio  (Ry)  required  for  calculating  Klqrg  from  K|_£, 

using   Equation  4.15.   (Note:   Critical  volume  for  ethylene   gas  =   120 
cm'/mol.) 

9.  Calculate  K[_qrg  for  the  organic  compound  from 

KL0RG  -  *LE/RV 


Example:  The  hydraulic  data  for  the  Speed  River,  survey  date  8/24/78, 
reach  BD,  given  in  Table  3.1,  have  been  considered  In  this  illustrative 
example.   The   organic    conpound    selected    1s   methylene    chloride 
(CH?  Clj).  The  computations  as  per  the  above  step-by-step  procedure 
are  as  follows: 

_> 

1.  For  methylene   chloride,    Henry's    law   constant  -  3  x  10 

ata-m'/mol  (see  Appendix  ), 

2.  The  hydraulic  data  for  the  reach  BO  of  the  Speed  River  are  as  follows: 

Q  •  2.43  m'/s;  l  -  2380  m;  S  =  0.0011;  B  -  38.6  ■; 

Z   •  0. 414  m;  U  -  13,392.0  m/day;  IS  »  (2380)(0.0Oll)  =  2.613  ■ 

3.  River  water  temperature  -  19.0*C 

o  -  998.4052  kg/m5;  y  =  88.73  kg/m-day  -  1.027  centipoise 
v  •  i»/p  ■  0.0889  mJ/day 

4.  The  Molar  volume  (Vg)  of  ethylene  gas  (C^Ha)  by  the  LeBas 
method  (see  Table  S.5  for  incremental  volumes)  1s: 

Ve  -  2(14.8)  *   4(3.7)  =•  44.4  cn'/mol 
2(C)     4(H) 

5.  From  Equation  3.4, 

13.26  x  10  "' 


-  1.3773  x  tO'1   cmVs 


(1.027)1-  "•(4*-*),'!" 

•  (1.3773  x  10"*)(cmI/s)(86400)(s/day)/(10*)(cmI/m2) 

•  1,19  x  10""  m'/day 

6.  From  Equation  4.9, 

K|_E  -[(2-7)  (7.3206  x  lO")0'1'  (Q.0889)*1 ' ' '  (1.19  x  lO"-)1-" 
(38.6)'-*  (0.414)-"-11  (^gZ)**"  (2.618)"-"] 
•  4.28  m/day 

7.  The  critical  volume  of  methylene  chloride  (CHjCl^),  calculated  by 
Lydersen's  method  outlined  1n  Section  5.5,  1s  as  follows: 

VC  -  40  ♦  SS  ♦  2(49)  -  193  cmVraol 
-  CH2  2(-Cl] 

The  crltictl  volume  of  ethylene  gas  Is  120  cm'/mol. 

8.  Fran  Equation  4.15, 

RV  •  ( 193/120)  "■,l*  ■  1.16 

9.  The  liquid  film  coefficient  of  methylene  chloride  1s 

K(_    «  (4.28/1.16)  -  3.69  m/day 

ORG 
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V  CHEMICAL  STRUCTURE  AND  VOLATILIZATION 

5.1  Relationship  Between  Henry's  Law  Constant 
and  Volatilization 

The  ability  of  an  organic  compound  to  volatilize  from  water  is  governed  by  two 
specific  parameters;   the  vapour  pressure  (Pv)  of  the  compound  and  the 
solubility  of  the  compound  (C$).  These  are  related  as  shown  in  Equation  5,1  to 
give  the  compound  specific  Henry's  law  constant  (H),  which  1s  used  as 

H  =  PV/CS  (5.1) 

a  measure  of  the  chemical's  ability  to  volatilize  from  water.   Pv  is  usually 
given  with  units  of  atm  and  C$  in  mol/m3  which  results  in  the  Henry's  law 
constant  (H)  having  units  of  atm  m3/mol.   When  H  is  calculated   using 
Equation  5.1  Pv  and  C$  used  should  be  for  the  same  temperature  and  physical 
state  of  the  chemical  and  the  data  should  be  for  the  pure  state.  The  relationship 
of  Henry's  law  constant  to  the  transfer  resistance  in  the  liquid  film  was  reviewed 
in  Chapter  II.   In  general,  a  compound  with  a  Henry's  law  constant  of  less  than 
10~7  atm  m3/mol  is  regarded  as  nonvolatile,  from  10~3  to 
10~7  as  moderately  volatile  and  greater  than  10~3  as  volatile. 

Since  Henry's  constants  are  not  readily  available  in  tabulated  form,  they  are 
usually  calculated  from  data  on  vapour  pressures  (Pv)  and  solubilities.  Vapour 
pressure  and  solubility  data  can  often  be  found  in  literature  sources  such  as  Weast 
(1980-81),  Perry  and  Chilton  (1963),  Dean  (1979),  Sax  (1979),  and  Verschueren 
(1977).  If  vapour  pressure  and  solubility  data  are  not  readily  available,  they  can 
be  estimated  from  the  basic  properties  of  the  compound  of  concern,  as  outlined  In 
the  following  sections. 

5.2  Estimation  of  Vapour  Pressures 

The  vapour  prpssures  for  organic  compounds  under  environmental  conditions  can  be 
estimated  from  the  normal  boiling  points  or  from  the  structure  alone  if  boiling 
point  data  are  not  available. 
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Estimation  of  Vapour  Pressure 

When  Boiling  Point  Data  are  Available: 

Two  methods  of  estimating  vapour  pressures  applicable  to  environmental  needs  have 
been  described  in  detail  by  Grain  in  Lyman  et  al  (1982).  These  are  based  on  the 
correlation  equation  of  Antoine  and  of  Watson  .  Both  of  these  are  derived  from 
Integration  of  the  Clausius-Clapeylon  equation  shown  below: 

d(S,n  Pv)    AHv 

= :  (5.2) 

dt     AZRT2 

where  AHV  is  the  heat  of  vaporization  for  the  compound;  AZ,  R  are 
defined  below;  and  T  is  the  temperature  in  degrees  Kelvin. 

Vapour  Pressure  from  the  Equation  of  Antoine: 

'The  procedure  based  on  the  Antoine  relationship  is  only  applicable  to  compounds 
that  are  in  the  liquid  or  vapour  phase  at  the  temperature  of  interest.  The 
equation  relating  vapour  pressure  (Pv)  and  the  normal  boiling  point  (Tb; 
at  one  standard  atmosphere  pressure)  is: 


£n  Py  = 


AHVD  (Tb-C2)2r   ,      ,   ,  (5-3) 


AZRTb2 


(Tb-C2)   (T-C2)J 


where   AZ  is  the  dimensionless  compressibility  factor  for 

the  gas  at  the  normal  boiling  point 'and  is  assumed 
to  be  equal  to  0.97  j 

AHvb  is  the  heat  of  vaporization  at  the  normal  boiling 

point  with  units  of  cal/mol  which  can  be  evaluated  using 
the  equation 

AHvb  =  KF  (8.75  +  R  In  Tb)  (5.4) 

the  gas  constant  R  is  equal  to  1.987  cal/deg  mol  and  <F  is  a 
constant  based  on  the  dipole  moment  of  the  compound.  Kp 
factors  are  available  from  a  number  of  sources  (Fishtine  i963; 
Lyman  et  al,  1982).  Typical  ranges  for  some  classes  of  organic   _ 
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compounds  are: 

1.  Hydrocarbons:  0.97  -  1.01 

2.  Halocarbons  (fluoro,  chloro,  bromo  and  iodo):  1.00  -  1.05 

3.  Carbonyl  compounds  (esters,  ketones  and  aldehydes): 
1.01  -  1.14 

4.  Nitrogen  compounds  (amines,  nitriles  and  nitros):  1.01  -  1.16 

5.  Thiols  and  Sulfides:  1.01  -  1.05 

6.  Monofunctional  alcohols:  1.20  -  1.31 

7.  Polyfunctional  alcohols  (glycols,  etc.):  1.33  -  1.38 

8.  Ethers:  1.01  -  1.08 

When  a  value  for  Kp  is  not  available,  a  default  value  of 
1.06  is  used.  C2  is  a  constant  whose  value  can  be  estimated 
from  (Thomson,  1959), 

C2  =  18  +  0.19  Tb  (5.5) 

Vapour  Pressure  from  the  Modified  Watson  Correlation: 

During  the  integration  of  the  Clausius-Clapeyron  Equation  (5.2)  if  the  ratio 
AHV/AZ  is  not  assumed  to  be  constant  with  varying  temperatures  but  is 
expressed  as  a  function  of  temperature,  the  resulting  expression  should  be  more 
exact.  Grain  in  Lyman  et  al  (1982)  reports  that  Equation  5.6,  a  modification  of 
the  Watson  correlation,  is  satisfactory  with  Tpb  =  T/Tb. 

AHV  =  AHvb  (3  -  2  Tpb)m  (5.6) 

and  m  is  a  constant  (see  below).  Equation  5.6  is  substituted  into  Equation  5.2 
which  is  then  integrated  by  parts.   Grain  states  that  sufficient  accuracy  is 
obtained  even  when  only  the  first  term  is  used.  Thus,  the  equation  relating  Pv 
with  AHvb  becomes: 


)m 
in  Pv  =       Jl       -^-   -2(1(3-  2Tpb)m~l  in  Tpb[    (5.7) 


AHvb  (     (3-2Tpb)' 
AZRTb I      Tpb 


The  following  values  for  the  constant  m  are  chosen  on  the  basis  of  the  physical 
state  of  the  compound  at  the  environmental  temperature: 
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1.  liquids  m  =  0.19 

2.  solids  -  if  Tpb  >  0.6,  m  =  0.36 

-  if  0.6  >  Tpb  >  0.5,  m  =0.8 

-  if  Tpb  <  0.5,  m  =  1.19 

R  is  equal  to  1.987  cal/deg  mol ,  AHvb/Tb  is  calculated 
using  Equation  5.3  and  AZ  =  0.97. 

The  value  of  Pv  obtained  by  this  method  has  units  of  atmospheres. 

Vapour  Pressure  from  Boiling  Points 
Obtained  at  Reduced  Pressures: 

For  many  organic  compounds  the  boiling  points  have  been  obtained  at  reduced 
pressures  for  experimental  convenience  or  to  avoid  decomposition  of  the  compound. 
As  a  result,  adaptions  of  Equations  5.3  and  5.5  are  needed  to  allow  the  estimation 
"of  Pv  for  these  compounds. 

Grain  has  approximated  the  ratio  of  the  heat  of  vapourization  to  temperature  as 


AHvi  _  KF  [8.75  +  R(JLn  Tx  -  In   Pi)] 
T, 


(5.8) 


When  Equation  5.8  is  used,  the  expression  for  vapour  pressure  based  on  Antoine's 
equation  becomes 


£n  P  ■  £n  P.  + 


AHvl  (T,  -  C2)^  r 


ZbRT 


r     1 


(TrC2)   (T-C2)J 


(5.9) 


The  expression  for  vapour  pressure,  based  on  the  modified  Watson  correlation, 
becomes 


£n 


P  «  *nP  ♦  ^ij1  "  l3"2  (T/Tl)im(ZiV  2m  t3_2  (T/Tl)3m_1  te(T/Tl>|     (5.1 


Note  that  Equation  5.9  is  only  valid  for  liquids  while  Equation  5.10  is  valid  for 
liquids  and  sol  ids. 


46 


I 


Estimation  of  Normal  Boiling  Point 
from  Chemical  Structure: 

The  normal  boiling  point  of  a  compound  is  defined  as  the  boiling  point  at  one 
standard  atmosphere.  Boiling  point  data  are  usually  available  in  the  literature- 
However,  for  those  compounds  where  data  is  not  readily  available,  an  estimation 
procedure  is  necessary. 

Most  methods  for  the  estimation  of  boiling  point  are  based  on  correlations  of 
parameters  that  are  affected  by  chemical  structure  with  known  boiling  points* 
Usually,  the  correlations  are  developed  using  relatively  simple  compounds 
containing  one  functional  group.  Chlorcbenzene  is  considered  to  have  one 
functional  group,  the  chlorine  substituent  on  the  benzene  ring.  It  is  usually 
assumed  that  multifunctional  compounds  can  be  treated  by  simply  adding  up  the 
effect  of  the  various  functional  groups.  For  example,  the  boiling  point  of 
dichlorobenzene  which  has  two  chlorine  functional  groups  would  be  estimated  by 
doubling  the  effect  of  having  only  one  chlorine  functional  group  as  is  the  case 
with  chlorobenzene. 

The  estimation  procedure  described  herein  is  that  of  Meissner  (1949).  The  method 
is  based  on  the  correlation  of  boiling  point  with  chemical  type.  The  equation  for 
the  boiling  point  in  degrees  Kelvin  (K)  is: 

637  (Rd]1'*7  +  B 
Tb  =  (5.11) 

[P] 

Rq  is  the  molar  refraction  and  can  be  estimated  using  the  atomic  refraction 
contribution  method  (see  belcw).  P  is  the  parachor  of  the  compound  which  may  be 
regarded  as  the  molecular  volume  of  a  liquid  having  a  surface  tension  of  unity 
(Quayle,  1953).  It  can  be  estimated  by  an  updated  version  of  Sugden's  method 
(Mumford  and  Phillips;  1928,  1929).  B  is  a  constant  whose  value  depends  on  the 
class  of  the  compound. 
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Estimation  of  Rq  and  P:  £< 

Rq  and  P  are  estimated  by  adding  the  contributions  appropriate  for  each  atom,     $ 
bond  type,  or  structural  unit.  The  contributions  for  each  structural  feature  are 
given  a  number  of  literature  sources  (Perry  and  Chilton,  1963;  Lyman  et  al  1982). 
The  method  is  best  described  using  chlorobenzene  as  an  example. 

Chlorobenzene  (CeHsCl)  has  the  structure: 


<0^ci 


The  molar  refraction  (Rq)  is  estimated  as: 

Rq  =  6  (contribution  for  a  carbon  atom) 

+5  (contribution  for  a  hydrogen  atom) 

+3  (contribution  for  a  double  bond) 

+1  (contribution  for  a  chlorine  atom) 

thus  RD  =  6  (2.418) 
-  +5  (1.100) 

+3  (1.733) 

+1  (5.967) 
RD  =  31.174 

Rq  could  be  more  efficiently  calculated  by  adding  the  contribution  for  a  phenyl 
group  (C6H5)  to  the  contribution  for  a  chlorine  atom  (CI).  Thus 

Rq  =  1  (contribution  for  a  phenyl  group) 
+1  (contribution  for  a  chlorine) 

RD  *  25.207 

+  5.967 

RD   31.174 


For  the  parachor  (P)  the  same  techniques  can  be  used,  eg: 

P  =  6  (  9.2) 

+5  (15.4) 

+3  (19.0) 

+1  (55.0) 

P  =  244.2 

or,  by  using  the  contribution  for  a  phenyl  group, 

P  =  1  (190.0) 
+   (55.0) 
P  =   245.0 

Note  that  in  the  calculation  of  P  using  a  value  of  190  for  the  phenyT  group  a 
slightly  different  value  of  P  is  obtained.  This  slight  variation  is  acceptable 
within  the  expected  error  of  the  method  (see  below).  For  B»  the  literature  (Fterry 
and  Chilton,  1973;  Lyman  et  al ,  1982)  gives  a  value  of  -2500  for  aromatic 
hydrocarbons.  This  value  of  B  has  been  chosen  as  there  is  no  value  for  Et  for 
chlorinated  aromatic  hydrocarbons.  To  calculate  a  value  of  %,  one  substitutes 
the  values  of  Rq,  P  and  B  into  Equation  5.11  to  obtain 

Tb  =  398  K  =  125  C 

As  a  comparison,  Weast  and  Astle  (1980-81)  reports  value  of  the  normal  bailing 
point  for  chlorobenzene  as  132°C.  The  calculated  value  of  E5*C  is  wittrin 
5%  of  the  reported  value.  Generally,  this  method  has  an  average  error  of  2%  v*fth  a 
maximum  error  of  about  8%. 

5.3   Estimation  of  the  Aqueous  Solubility  of  a  Compound 

As  indicated  by  Equation  5.1,  Henry's  law  constant  is  inversely  proportional  to  the 
aqueous  solubility.  Solubilities  of  organic  compounds  vary  aver  a  wide  range;  with 
most  organic  compounds  having  solubilities  from  1  to  100,00©  ug/k  although  soraze 
organics  have  solubilities  as  low  as  0.001  ug/L  while  others  are  infinitely 
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soluble  1n  water.  Due  to  this  extremely  wide  range  of  solubilities  the  value  of 
the  Henry's  constant  -  which  governs  the  volatilization  rate  -  is  highly  dependent 
on  the  solubility. 

Literature  values  for  the  solubility  of  organic  compounds  are  usually  readily 
available.  However,  a  method  of  estimating  the  solubility  of  organics  in  water  is 
needed  for  compounds  for  which  data  are  not  obtainable  or  for  which  none  has  been 
determined. 

Solubility  From  Octanol-Water  Partition  Coefficients: 

A  number  of  correlation  relationships  between   the   octanol-water   partition 
coefficient  (Kow)  and  the  aqueous  solubility  (C$)  have  been  presented  in  a 
recent  publication  Lyman  et  al,  (1982).  Table  5.1  presents  a  selected  number  of 
these  correlations.   An  example  of  the  use  of  one  of  these  correlation  equations 
will  be  presented  during  the  discussion  of  the  octanol-water  partition  coefficient. 

Octanol-Water  Partition  Coefficient 
from  Chemical  Structure: 

Lyman  et  al  (1982)  present  a  review  of  the  available  methods  for  the  estimation  of 
Kow.  The  methods  include  those  based  on  fragment  constants  which  are  summed  to 
give  the  octanol-water  partition  coefficient,  methods  based  on  the  estimation  of 
activity  coefficients,  by  calculation  of  free  energies  of  solvation,  via 
topological  indices,  and  by  linear  regressions  on  to  other  physical-chemical 
parameters.  The  method  which  best  meets  the  criterion  of  accuracy,  ease  of 
calculation,  and  applicability  to  a  wide  range  of  compounds  is  that  of  Hansch  and 
Leo  (1979).  This  method  makes  use  of  fragment  constants  (f)  and. structural  factors 
(F)  to  estimate  Kow.  Values  of  f  and  F  can  be  found  in  a  number  of  sources 
(Hansch  and  Leo,  1979;  Lyman  et  al ,  1982).  In  addition,  Hansch  and  Leo  (1979) 
provide  a  listing  of  Kow  values  for  thousands  of  compounds. 

Calculations  of  the  octanol-water  partition  coefficient  are  carried  out  in  terms  of 
log  Kow  using  the  equation 

Log  Kow  -  sum  of  fragments  (f)  +  factors  (F)  (5.12) 
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Table  5.1 

EQUATIONS  RELATING  AQUEOUS  SOLUBILITY 
TO  OCTANOL -WATER  PARTITION  COEFFICIENT 


tn 


Equation 


log  S  -  -1.12  log  K   +  7.30  -  0.015  t+ 

ow  m 


log  1/S  - -1.339  log  K   -  0.978 

OW 


Ref 


log  S  -  -2.38  log  K       +  12.90 

b 

log  S  =  -0.9874  log  K       -  0.0095  t     +  0.7178 

ow                     m 

b 

log  S  =  -0.88  log  K       -  0.01  t     -  0.012 

ow                 m 

b 

log  S  =  -0.962  log  K      +  6.50 

ow 

b 

Units  of 

Solubility  (s) 

umol/L 

mol/L 

UlTOl/L 

umol/L 

mol/L 

7  u  mol/L 

Classes  of  Compounds 

All  types;  especially  aromatlcs 
and  chlorinated  hydrocarbons 

Alcohols,  Ketones,  Esters,  Ethers, 
Alkyl  halids,  Alkynes,  althenes, 
Aromatics,  Alkanles 

Phosphate  esters 

Halogenated  benzene  compounds 

Polycyclic  Aromatic  Hydrocarbons 

Halogenated  organics  containing  1  or 
2  carbons  only 


Banarjete  et  al  (1980) 

Lyman  et  al  (1982)  and  references  cited 

tm  is  the  melting  point  of  the  compound.  If  tm  is  <25°C,  a  default  value  of  25°C  1s  used. 


m 


Although  a  complete  description  of  the  method  is  beyond  the  scope  of  this  repc., 
an  example  is  provided  in  order  to  illustrate  the  method.  The  following  example  -s 
taken  from  Hansch  and  Leo  (1979).  The  structure  of  lindane  (CgHeCle) 
1s: 


The  fragments  (f)  necessary  for  calculation  of  Kow  are: 

6  (0.20  for  each  carbon); 
6  (0.23  for  each  hydrogen); 
and  6  (0.06  for  each  chlorine) 

The  Factors  (F)  necessary  for  the  calculation  of  Kow  are: 

(6-1)  (-0.09  for  bonds  in  a  ring); 
and  (6-1)  (0.28  for  vicinal  chlorines). 

By  substituting  these  values  in  Equation  5.12,  the  calculated  value  of  log  Kow 
for  lindane  is  3.89.   The  observed  value  reported  by  Hansch  and  Leo  (1979)  was 
3.72.  When  the  calculated  log  Kow  is  substituted  into  the  first  equation  of 
Table  5.1,  an  aqueous  solubility  of  0.018  mol/m3  is  obtained.  This  compares 
reasonably  well  with  the  value  of  0.034  mol/m3  reported  by  Verschueren  (1977). 

5.4   Henry's  Constants  for  Compounds  of  Concern  in  Ontario 

The  available  data  for  the  compounds  of  concern  to  the  Ministry  are  listed  in 
appendix.  In  this  appendix,  the  compounds  have  been  divided  into  three  groups: 


High  volatility:  H  >  10" 3  atm  m3/mol; 

Moderate  volatility:  10"3  atm  m3/mol  >  H  >  10"7  atm  m3/mol); 

Low  volatility:  H  <  10" 7  atm  m3/mol). 


and 
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The  values  of  the  Henry's  law  constant  which  were  used  to  make  the  classification 
are  literature  values  from  the  references  cited  or  estimated  using  literature  data 
on  boiling  points,  vapour  pressures,  solubilities,  etc.  Where  data  was 
unavailable,  an  estimation  of  the  Henry's  law  constant  was  made  using  the 
techniques  of  the  preceding  sections  of  this  chapter.  Within  each  group,  the 
compounds  are  classified  into  chemical  families. 


5.5   Estimation  of  Molar  and  Critical  Volumes 
5.5.1  Molar  Volume  from  Chemical  Structure: 

The  molar  volume  (V3)  can  be  calculated  from  a  knowledge  cf  structure  alone,  by 
the  LeBas  method.  The  method  has  been  described  in  texts  such  as  Perry  and  Chilton 
(1973),  Reid  et  al,  1977  and  Lyman  et  al  (1982).  In  the  LeBas  method,  Vg  is 
calculated  by  summing  the  additive  volume  increments  shown  in  Table  5.2  for  each 
atom  and  structural  feature  in  the  molecule.  The  procedure  is  best  described  using 
chlorobenzene  (C6H5CI)  as  an  example.  The  structure  cf  chlorobenzene  is 


C^ci 


From  an  examination  of  this  structure  the  LeBas  molar  volume  is  estimated  in  the 
following  manner: 

Vg  =    1  (increment  for  chlorine) 

+   6  (Increment  for  carbon 

+   5  (increment  for  hydrogen) 

+   1  (increment  for  6-membered  ring) 
and 

VB  =  24.6 

+   6   (14.8) 

+   5   (3.7) 

+       1       (-15.0) 

V*b  ~  116.9  cmVrnol 
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5.5.2  Critical  Volume  from  Chemical  Structure: 

Critical  volume  is  defined  as  the  volume  of  a  mole  of  a  compound  at  the  highest 
temperature  that  the  compound  can  exist  as  a  liquid.  Lydersen  has  developed  a 
method  of  calculating  tfoe  critical  volume  (Vrj)  from  only  a  knowledge  of  the 
structure  of  the  compound.  Lydersen 's  method  has  been  presented  in  texts  such  as 
Reid,  et  al  (1977)  and  Lyman  et  al  (1982).  The  method  is  based  on  the  summation  of 
volume  increments  for  esach  of  the  structural  features  of  the  organic  compound.  A 
listing  of  the  structural  features  and  the  associated  volume  increments  (AV)  is 
given  in  Table  5.3.  Ths  equation  for  Vc  is 

Vq  =  40  +  £AV 

where  the  critical  velume  has  the  units  of  ctn'/mol.  The  procedure  for 
calculating  Vr,  is  best  described  using  an  example,  specifically  chlorobenzene, 
which  has  the  structure  shown  in  section  5.5.1. 


From  an  examination  off  the  structure  in  conjunction  with  the  data  given  in 
Table  5.3, 

Vc  ■  40 

+  AV  (chlorine) 

+  5  AV  (HC=  in  a  ring) 

+  AV  (C=  in  a  ring) 

VC  ■   40 
49 
5(37) 
♦  36 
Vc  ■  310  cnrVmol 
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Table  5.2 


ADDITIVE  VOLUME  INCREMENTS  FOR  CALCULATING 
THE  LEBAS  MOLAR  VOLUME 


Atom 

Increment 

Atom 

Increment 

{emVmol) 

(cm'/mol) 

C 

14.8 

Br 

27.0 

H 

3.7 

CI 

24.6 

0  (except  as  noted  below) 

7.4 

F 

8.7 

In  methyl  esters  and  ethers 

9.1 

I 

37.0 

In  ethyl  esters  and  ethers 

9.9 

s 

25.6 

In  higher  esters  and  ethers 

11.0 

In  acids 

Joined  to  S,  P,  N 
N 

Double  bonded 
In  primary  amines 
In  secondary  amines 

12.0 
8.3 

Ring 
3-Membered 
4-Membered 

-  6.0 

-  8.5 

15.6 
10.5 
12.0 

5-Membered 
6-Membered 
Naphthalene 
Anthracene 

-11.5 
-15.0 

-30.0 
-47.5 

Source:  Reid,  et  al  (1977) 


55 


Table  5.3 


LYDERSEN'S  INCREMENTS  FOR  CALCULATING  CRITICAL  VOLUMES 


a,b,c 


Increment 

AV 

Nonring  increments 

-CH, 

55 

1 

-CH3 

55 

! 
-CH 

1 

51 

1 

-C- 

| 

41 

=CH2 

45 

1 
=CH 

45 

1 

=C- 

36 

=C= 

36 

sCH 

<36)c 

=C- 

(36) 

Ring  increments 

-CH,- 

44.5 

I 

-CH 

1 

46 

X 

-     (31) 

Increment 


AV 


Halogen  increments 

-F 

-CI 

-Br 

-I 
Oxygen  increments 

-OH  (alcohols) 

—OH  (phenols) 

—0—  (nonring) 

— O—  (ring) 

I 
— C=0  (nonring) 

I 
-C=0  (ring) 

i 
HC=0  (aldehyde) 

-COOH  (acid) 

-COO-  (ester) 

=0  (exc.  for  above 
combinations) 

Nitrogen  increments 

-NH2 

I 
— NH  (nonring) 


IB 

49 

(70) 

95 

(18) 
(3) 

20 
(8) 

60 

(50) 

73 
80 
80 

11 


28 


=CH 


37 


-C- 

36 

-C- 

36 

Ifur  increments 

-SH 

55 

— S—  (nonring) 

55 

— S—  (ring) 

(45) 

=S 

(47) 

-NH(ring) 

I 

— N—  (nonring) 

I 
-N-  (ring) 

-CN 

-NOj 


(27) 
(42) 
(32) 


(78) 


*.     Th«rt  are  no   Increments  for  hydrogen. 

b.  Alt  bonds  shown  as  free  *re  connected  to  atoms  other  than  H. 

c.  Values  In  parentheses  are  based  upon  too  few  experimental 
values  to  be  rel  (able  . 


Source:     Re1d.  et  al   (1977). 
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Vt         CONCLUSIONS  AND  RECOMMENDATIONS 

Volatilization  of  organic  compounds  in  streams  and  rivers  has  been  investigated 
through  an  analysis  of  data  on  ethylene  and  dye  tracer  and  development  of 
predictive  relationships.  The  factors  affecting  volatilization,  and  conceptual 
models  of  transport  at  the  air-water  interface  have  been  presented.  The 
relationship  of  liquid  film  coefficients  to  bulk  flow  channel  parameters  and 
characteristic  properties  of  organic  compounds,  reported  in  the  literature,  have 
been  reviewed.  The  liquid  film  coefficients  of  ethylene  for  provincial  streams  and 
rivers  have  been  determined  by  making  use  of  data  on  ethylene  and  rhodamine  WT  dye 
tracers  which  were  collected  in  the  provincial  streams  and  rivers.  The  development 
of  generalized  relationships  among  the  liquid  film  coefficients  and  bulk  flew 
channel  characteristics  has  been  accomplished  through  statistical  and  dimensional 
analysis  techniques. 

The  predictive  capabilities  of  the  relationships  have  been  evaluated  by   a 
comparative  study  of  the  observed  and  the  predicted  values  of  K[_  for  ethylene 
gas  applicable  to  the  provincial  streams  as  well  as  other  streams  for  which  data 
are  reported  in  the  literature  (Rathbun  and  Tai,  1981;  1982).  A  sensitivity 
analysis  of  various  parameters  appearing  in  the  relationships  has  been  presented. 

The  relation  of  liquid  film  coefficient  (K|_)  of  an  organic  compound  to  the 
molecular  weight  as  well  as  critical  volume  (or  molecular  diameter)  of  the  compound 
have  been  evaluated  with  the  aid  of  laboratory  experimental  data  reported  in  the 
literature.  From  this  evaluation,  relationships  for  computing  K[_  values  of 
organic  compounds  from  the  K^  values  for  ethylene  gas,  have  been  presented. 

The  computation  of  liquid  film  coefficient  of  an  organic  compound  using  the 
relationships  developed  herein  is  outlined  in  a  step-by-step  procedure,  and 
illustrated  by  an  example. 

The  relation  between  chemical  structure  and  volatilization  has  been  assessed  for  a 
number  of  organic  compounds  of  concern  in  the  province.  The  compounds  have  been 
grouped  into  a  number  of  classes  which  are  expected  to  behave  similarly  in  streams 
and  rivers  with  respect  to  volatilization  based  on  the  Henry's  law  constants. 
Methods  of  computing  Henry's  law  constants,  molar  volumes  and  critical  volumes  of 
organic  compounds  have  been  outlined. 
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The  major  conclusions  of  this  investigation  are  as  follows: 

1.  The  liquid  film  coefficients  of  ethylene  in  the  provincial  streams  and 
rivers  considered  in  this  study  are  in  the  range  0.07  to  8.07  m/day  at  the 
observed  average  temperatures  (or  1.07  -  7.69  m/day  at  20°C). 

2.  The  predictive  abilities  of  the  relationships  among  liquid  film  coefficient 
and  bulk  flow  characteristics  derived  from  dimensional  analysis  procedures 
are  better  than  the  statistical  correlation  relationships. 

3.  The  relation  obtained  by  Rathbun  and  Tai  (1982)  underestimates  the  liquid 
film  coefficients  of  ethylene  at  20°C  for  the  provincial  streams  and 
rivers. 

4.  Application  of  Equation  4.10  developed  herein  from  dimensional  analysis 
procedure  to  other  streams  in  the  literature  {Rathbun  and  Tai,  1981) 
results  in  an  overestimation  of  the  liquid  film  coefficients  for  ethylene 
at  20°C.  However,  the  coefficients  were  indirectly  calculated  (but  not 
directly  measured). 

5.  Equations  4.9  and  4.10  can  be  utilized  to  estimate  liquid  film  coefficients 
for  ethylene  gas  in  provincial  streams  and  rivers.  These  relationships  do 
not  include  wind  velocity  and  hence,  the  computed  values  may  involve 
approximations.  These  limitations  should  be  given  due  consideration  in  the 
interpretation  and  application  of  the  results. 

6.  The  ratios  (<l  for  ethylene  gas/K|_  for  organic  compound)  obtained 

from  experimental  data  are  in  reasonable  agreement  with  the  values  computed 
by  using  the  critical  volumes  of  the  compounds  (Equation  4.15)  but  differed 
from  those  obtained  by  using  molecular  weights  (Equation  4.16)  for  organic 
compounds  having  Henry's  law  constants  >  10   atm  -  m  /mol.   Thus, 
the  ratios  obtained  from  Equation  4.15  are  applicable  to  calculate  Kj_ 
values  of  highly  volatile  organic  compounds  from  those  of  ethylene  gas. 
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As  outlined  above,  an  evaluation  of  the  predictive  capabilities  of  the 
relationships  among  K|_,  bulk  flow  channel  parameters  and  properties  of  chemical 
compounds  indicates  that  the  relationships  may  not  be  universally  applicable  to 
streams  and  rivers  of  differing  hydraulic  and  water  quality  characteristics,  as 
well  as  for  different  organic  chemicals  of  public  health  concern.  In  general, 
there  is  a  deficiency  in  the  data  required  for  development  and  validation  of 
generalized  relationships  among  Kj_,  channel  hydraulic  parameters  and  properties 
of  chemicals.  The  effect  of  wind  velocity  is  not  included  in  the  relationships 
because  of  a  lack  of  field  data.  Based  on  the  knowledge  obtained  in  this  study, 
the  following  recommendations  for  futher  work  are  made: 


1.  The  development  of  modified  relationships  for  K|_  to  include  wind 
velocity  through  the  dimensional  analysis  procedure,  is  recommended. 

2.  It  is  recommended  to  test  the  applicability  of  the  relationships  to 
chemical  compounds  other  than  ethylene  as  data  become  available. 

3.  Detailed  field  data  collection,  including  the  channel  hydraulic  parameters 
(discharge,  width,  depth,  velocity,  bed  slope)  and  wind  velocity 
measurement,  are  recommended  to  be  undertaken  in  future  studies  on 
volatilization  of  organic  compounds  of  public  health  concern  from  natural 
water  bodies. 

4.  Investigation  of  relationships  between  K[_  and  molecular  properties  for 

chemical   compounds   with   Henry's    constants     lower         than 

-3       3 

10  "  atm  -  m  /mol  (i.e.,  moderate  and  low  volatile  compounds)  is 

recommended. 
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APPENDIX 


COMPOUNDS  OF  HIGH  VOLATILITY 


Ha 
atm  m/rnol 


atiYi 


mol/nr 


Other  Data 


Hydrocarbons 
Benzene 
Butyl ene 
Cyclohexane 
Ethylene 
Ethyl  benzene 

Isopropylbenzene 
n-Hexane 
Natural  gas 
Propylene 
Styrene  (monomer) 
Toluene 


4.4  x  10"J  (1) 

>  10"3  (7) 

0.15 

>  8.6 
8.7  x  10"3 

5.7  x  10"3  (6) 

9.8  x  10"3  (5) 

1.4 
»  10"3  (7) 


-3 


2.1  x  10 


4.3  x  10"J  (5) 

5.2  x  10"3  (1) 

5.7  x  10"3  (6) 

6.6  x  ID"3  (4) 


0.10  (5) 
>  40  (4) 
1.25  x  10"2  (4) 


4.2  x  10'°  (5) 
0.16  (5) 

10  (12) 
6.6  x  10"3  (9) 


0.65  (5) 
4.7  (4) 
1.43  (4) 

0.42  (5) 
0.11  (5) 

4.75  (b) 
1.54  (9) 


3.7  x  10"^  (4) 


64 


Xylene 
ortho 
meta 
para 


5.3  x  10"J  P  25  C 
4.2  x  10'3  (7) 
4.6  x  10"3 


Polycylic  Aromatic  Hydrocarbons  (PAHs) 


Anthracene 

Naphthalene 

Oxygenated  Hydrocarbons 

Pesticides 

Aldrin 

alpha-Hexachlorocyclohexane 

Toxaphene 

Halogenated  Hydrocarbons 

BromoiTie  thane 
(methyl  bromide) 

Carbon 
Tetrachloride 

Chlorobenzene 
Chloroethane 


1.4  x  10~J  (6) 
1.18  x  10"3  (4) 


2.1  x  10""3  (6) 
2.0  x  10"3  (6) 
6.3  x  10~2  (6) 


0,25 

£,5  *  10^8  (6) 

2.6  x  10r3  (1) 

4:Q  *   10-3  ($) 

1.5  x  10-2  (6) 


8.7  x  10"J  9   25  C  (3)    1.64  *  25  C  (3) 


7.9  x  10"J  (5) 
8.6  x  10"3  (5) 


M  (I) 


1.87  (5) 


9.5  (5) 
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Chloroform 

Chloromethane 

3-Chloropropene 

1,2-Di  chlorobenzene 

1,3-Dichlorobenzene 
(m-di chlorobenzene) 

1,4-Di chlorobenzene 
(p-di chlorobenzene) 

Dichlorodi  f luoromethane 

1,1-Dichloroethar.e 

1,1-Dichlorethene 

(vinylidene  chloride) 

cis-l,2-Dichloroethene 

trans-l,2-Dichloroethene 

1,2-Di chloropropane 

Hexachloroethane 

Methylene  chloride 

PCBs 
Aroclor  1248 

Aroclor  1254 

Aroclor   1260 

Pentachloroethane 


3.4     x  10"3  (6) 

31 
1.1     x     ltf2  (10) 
1.7     x  10"3  (6) 
2.7     x   10"3  (6) 


2.1     x  10";  (6) 
3.4     x   10"J  (9) 

2.4 

5.1     x   10"3  (6) 

1.7     x   10'1   (6) 


7.5  x   10       (10) 

6.5  x  10"3  (10) 

2.1  x  10"3  (1) 

1.1  x   10"3  (6) 

3.0  x   10"3  (4) 


3.51  x  10"J   (3) 

2.76  x  10"3  (3) 

7.13  x  10-3  (3) 

•2.4     x  10-3  (10) 


5.0     (5) 


0.16  (5) 


5.6     (5) 


2.3     (5) 


0.46   (4) 


155         (4) 


Perchloroethylene 
(tetrachloroethene) 

1,1,1,2-letrachloroethane 

1,1,1-Trichloroethane 

Trichloroethene 

Trichlorof luoromethane 

2,3  Dichloropropene 

cis-1 ,3-Dichloropropene 

trans-1 ,3-Dichloropropene 

Vinyl  chloride 
(chloroethene) 


8.3     x   10"°  (4) 


2.6     x  1Q~J  (10) 

1.8     x   10"^  (4) 
3.6     x  1D~J  (6) 

1.0     x   10"2   (6) 

0.11   (6) 

3.6  x   10"3   (10) 
2.3     x   10'3  (10) 

1.7  x   10"3  (10) 

2.4     (4) 
6.4     (6) 


2  x   10~2  (4) 


2.4     (4) 


3.4     (4) 


1.44  (4) 


Nitrogen  Compounds 
Methylamine 


8.3  x  10 


3.1  (5) 


37.2  (5) 


Miscellaneous  Compounds 
Carbon  disulfide 


1.1  x  10" 


0.34  (5) 


30.2  (b) 


C0MP00NDS  OF  MODERATE  VOLATILITY 


Hydrocarbons 
Acetylene 


7.5  x  lO-^  (7) 


40    (5) 
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1,3-Butadiene 


2.4     x  10 


-4 


3.3     x   10'J   (5) 


13.6  (5) 


Polycycllc  Aromatic  Hydrocarbons   (PAHs) 


Benzo(a)pyrene 
Acenaphthene 
Fluorene 
Phenanthrene 

Oxygenated  Hydrocarbons 

Acetone 

Acrolein 

n-Butanol 

n-Butyraldehyde 

Cresols 

ortho 

meta 

para 
Dimethyl  ether 
Dimethyl  phthalate 
4,6-Dinitro-o-cresol 
Ethylene  Oxide 


1.3  x  10'°  (11) 

1.9  x  10'4  (6) 

2.1  x  10"4  (6) 

1.3  x   10"4  (6) 


1.5     x  10"3  (7) 
9.7     x  10"5  (6) 


-6 


5.6     x   10 


1.0  to  2.2  x  10       (7) 


1.1  x  10~6  (7) 
2.4  x   10'7  (7) 

8.6  x   10"6  (11) 
3.3  x   10"4  (7) 

4.2  x   10"7  (6) 

1.7  x   1Q~7  (6) 
7.7  x   10"b  (7) 


5.8     x  10"3  (5) 


5.3     x   ID"3  (5) 


5.0  (12) 


1040  (5) 
430  to  985  (5) 


1.4   (5) 


2-Ethylhexanol 

Isophorone 

Isopropanol 

Methyl   Ethyl  Ketone 

Methyl    Isobutyl  Ketone 

2-Nitrophenol 

Pentachlorophenol 

Phenol 

Vinyl  Acetate 


2.1 

x  10  "  (7) 

5.8 

x  10"6 

9.3 

x  10"6  (7) 

2.0 

-b 
x   10  D 

4.7 

-5 
X   10  J 

7.6  x  10"3  (6) 
3.4  x   10"6  (4) 

2.7  x  10"7  (6) 

3.8  x   10~4 


3.1   (5) 
5.0     x  10"4  (b) 
4.2     x   10"2  (5) 

0.10  (5) 
7.9     x   10'3  (5) 


1.8     x  10~'   (4) 


1.1     x   10"1   (5) 


87   (5) 

4900  (b) 
170  (5) 


5.3  x  10"2  (4) 


290  (5) 


Pesticides 

Chlordane 

DDT 

(4,4'-DDl) 

Dieldrin 

Lindane 

Methyl  parathion 

Mi  rex 


3.7     x  10"H  (7) 


3.4  x  10'5  (6) 

1.7  x   10"7  (6) 
3.2  x   10"7  (6) 

7.8  x   10"5  (11) 

1.9  x   10"4  (11) 


Halogenated  Hydrocarbons 
Bis(2-chloroethyl )ether 
Bis(2-chloroisopropyl )ether 


-5 


1.3     x   10 

1.1     x  10~4   (6) 


9.3     x  10"*  (5) 


71.3  (5) 
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1,2-01 chloroethane 

Ethylene  Dlbromide 

PCBs 
Aroclor  1242 

Tribromo methane 

1,2,3-lMdhloropropM! 

1,1,2,2-Tetrachloroethane 
1,1, 2-1 ri chloroethane 

Nitrogen  Compounds 

Acrylonitrile 

Ammonia 

Carbazole 

Nitrobenzene 

Hydrazine 

Miscellaneous  Compounds 

Benzo(b)thiophene 

Dibenzothiophene 


9.6     x  10"4  (10) 
6.6     x  10"4  (4) 


4.9  x  10~H  (6) 

5.9  x   10'4 

3.1  X  10"4  (10) 

4.2  x  ltf4  (6) 
7.8  x  10'4  (6) 


6.3  x   10"J  (6) 

2.7  x   10"4 

1.6  x   10"4  (11) 

1.1  x  10"b  (6) 

1.8  x   10"5  (7) 


2.7     x   10"4  (11) 
4.4     x   10"4  (11) 


7.4     x  10""3  (5) 


8.7     (5) 


2.1     x  10"'  (5) 


12.6  (5) 


32000  (5) 


COMPOUNDS  OF  LOW  VOLATILITY 


Poly cyclic  Aromatic 

Hydroca 

rbons 

(PAHs) 

Benz(a)anthracene 

1.3     x  10"7  (11) 

Oxygenated  Hydrocarbons 

Ethanol 

_7 

<   10   ' 

5.8 

x   10'2  (5) 

Ethanolamine 

9.1     x   10"8   (7) 

5.3 

x   10"4  (5) 

Ethylene  glycol 

«   10"7 

6.6 

x   10"5 

mlsdble 

Methanol 

«   10"7 

1.2 

x   10"1 

mlscible 

Nitrogen  Compounds 

Benzoff  )qu1noline 

9.5     x   10"9  (11) 

E-Caprolactain 

2.6     x   10"10   (7) 

1.3 

x   10'6  (5) 

Di  benzo (c , g)carbazol 

e 

4.7     x   10"9  (11) 

Morpholine 

4     x   10"8  (7) 

1.1 

x   10"2  (5) 

Urea 

<   10'7   (7) 

Nitroglycerine 

4.1     x   10"8  (7) 
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Qulnollne  2-5  *  10~  (») 


a:  For  20°C  unless  otherwise  noted 

1.  Mackay  and  Yuen  (1983) 

2.  Banerjee  et  al  (1980) 

3.  Mackay  and  Leinonen  (1975) 

4.  Mackay  et  al  (1979) 

5.  Verchueren   (1977) 

6.  McCarty   (1980)  calculated  values 

?.  P     and/or  C     estimated  using  the  methods  of  Chapter  b 

8.  sSx   (1979)   S 

9.  Arbuckle   (1983) 

10.  Dilling  (1977);  some  values  were  estimated.  See  reference  for  details 

11.  Smith  and  Bomberger  (1978) 

12.  Weast  and  Astle  (1980-81) 
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